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Abstract 



In the last 20 years, loop quantum gravity, a background independent approach to unify general 
relativity and quantum mechanics, has been widely investigated. The aim of loop quantum 
gravity is to construct a mathematically rigorous, background independent, nonperturbative 
quantum theory for the Lorentzian gravitational field on a four-dimensional manifold. In this 
approach, the principles of quantum mechanics are combined with those of general relativity 
naturally. Such a combination provides us a picture of "quantum Riemannian geometry", which 
is discrete at a fundamental scale. In the investigation of quantum dynamics, the classical 
expressions of constraints are quantized as operators. The quantum evolution is contained in 
the solutions of the quantum constraint equations. On the other hand, the semi-classical analysis 
has to be carried out in order to test the semiclassical limit of the quantum dynamics. 

In this thesis, the structure of the dynamical theory in loop quantum gravity is presented ped- 
agogically. The outline is as follows: first we review the classical formalism of general relativity 
as a dynamical theory of connections. Then the kinematical Ashtekar-Isham-Lewandowski rep- 
resentation is introduced as a foundation of loop quantum gravity. We discuss the construction 
of a Hamiltonian constraint operator and the master constraint programme, for both the cases of 
pure gravity and matter field coupling. Finally, some strategies are discussed concerning testing 
the semiclassical limit of the quantum dynamics. 
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1 Introduction 



1.1 Motivation of Quantum Gravity 

The current view of physics is that there exist four fundamental interactions: strong interaction, 
weak interaction, electromagnetic interaction and gravitational interaction. The description for 
the former three kinds of forces is quantized in the well-known standard model. The interactions 
are transmitted via the exchange of particles. However, the last kind of interaction, gravitational 
interaction, is described by Einstein's theory of general relativity, which is a classical theory 
which describes the gravitational field as a smooth metric tensor field on a manifold, i.e., a 4- 
dimensional spacetime geometry. There is no h and hence no quantum structure of spacetime. 
Thus there is a fundamental inconsistency in our current description of the whole physical 
world. Physicists widely accept the assumption that our world is quantized at fundamental level. 
So all interactions should be brought into the framework of quantum mechanics fundamentally. 
As a result, the gravitational field should also have "quantum structure". 

Throughout the last century, our understanding of nature has considerably improved from 
macroscale to microscale, including the phenomena at molecule, atom, sub-atom, and elemen- 
tary particle scale. The standard model of particle physics agrees with all present experimental 
tests in laboratory (see e.g. I1158II ). However, because unimaginably large amount of energy 
would be needed, no experimental tests exist for processes that happen near the Planck scale 
£p = (Gn/cy^ ~ 10~^^cm and tp = (GH/cy^ ~ lO"'^-^^, which are viewed as the most 
fundamental scales. The Planck scale arises naturally in attempts to formulate a quantum the- 
ory of gravity, since £p and tp are unique combinations of speed of light c, Planck constant fi, 
and gravitational constant G, which have the dimensions of length and time respectively. The 
dimensional arguments suggest that at Planck scale the smooth structure of spacetime should 
break down, and therefore the well-known quantum field theory is invalid since it depends on 
a fixed smooth background spacetime. Hence we believe that physicists should go beyond the 
successful standard model to explore the new physics near Planck scale, which is, perhaps, a 
quantum field theory without a background spacetime, and this quantum field theory should 
include the quantum theory of gravity. Moreover, current theoretical physics is thirsting for a 
quantum theory of gravity to solve at least the following fundamental difficulties. 

• Classical Gravity - Quantum Matter Inconsistency 

The equation relating matter and the gravitational field is the famous Einstein field equa- 
tion: ^ 

Raflig] - :^R[g]ga/i = KT^fiigl (1) 

where the left hand side of the equation concerns spacetime geometry which has classical 
smooth structure, while the right hand side concerns also matter field which is funda- 
mentally quantum mechanical in standard model. In quantum field theory the energy- 
momentum tensor of matter field should be an operator- valued tensor T^^. One possible 
way to keep classical geometry consistent with quantum matter is to replace T^^lg] by 
the expectation value < T^^lg] > with respect to some quantum state of the matter on a 
fixed spacetime. However, in the solution of this equation the background gajs has to be 
changed due to the non- vanishing of < fafj[g] >. So one has to feed back the new metric 
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into the definition of the vacuum expectation value etc. The result of the iterations does 
not converge in general fTOl . On the other hand, some other arguments show that such a 
semiclassical treatment may violate the principle of superposition in quantum mechanics 
ll55l . This inconsistency motivates us to quantize the background geometry to arrive at an 
operator formula also on the left hand side of Eq.©- 

• Singularities in General Relativity 

Einstein's theory of General Relativity is considered as one of the most elegant theories 
in the 20th century. Many experimental tests confirm the theory in the classical domain 
ni59l| . However, Penrose and Hawking proved that singularities are inevitable in general 
spacetimes with matter satisfying certain conditions in, by now well known, singularity 
theorems (for a summary, see [ f9T| [ ll56ll ). Thus general relativity as a classical theory 
breaks down in certain regions of spacetime in a generic way. One naturally expects that, 
in extra strong gravitational field domains near the singularities, the gravitational theory 
would probably be replaced by an unknown quantum theory of gravity. 

• Infinities in Quantum Field Theory 

It is well known that there are infinity problems in quantum field theory in Minkowski 
spacetime. In curved spacetime, the problem of divergences is even more complicated, 
since the renormalization process in curved spacetime is ambiguous, the expectation value 
of stress tensor can be fixed up to some local curvature terms, and it also depends on a 
fundamental scale of spacetime. Although much progress on the renormalization have 
been made |[921[|157ll , a fundamentally satisfactory theory is still far from reaching. So it 
is expected that some quantum gravity theory, playing a fundamental role at Planck scale, 
could provide a natural cut-off to cure the infinities in quantum field theory. The situation 
of quantum field theory on a fixed spacetime looks just like that of quantum mechanics 
for particles in electromagnetic field before the establishing of quantum electrodynamics, 
where the particle mechanics (actress) is quantized but the background electromagnetic 
field (stage) is classical. The history suggests that such a semi-classical situation is only 
an approximation which should be replaced by a much more fundamental and satisfactory 
theory. 

1.2 Purpose of Loop Quantum Gravity 

The research on quantum gravity is quite active. Many quantization programmes for gravity 
are being carried out (for a summary see e.g. [146J). In these different kinds of approaches. 
Among these different kinds of approaches, the idea of loop quantum gravity finds its roots in 
researchers from the general relativity community. It follows closely the motivations of general 
relativity, and hence it is a quantum theory born with background independence. Roughly 
speaking, loop quantum gravity is an attempt to construct a mathematically rigorous, non- 
perturbative, background independent quantum theory of four-dimensional, Lorentzian general 
relativity plus all known matter in the continuum. The project of loop quantum gravity inherits 
the basic idea of Einstein that gravity is fundamentally spacetime geometry. Here one believes 
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in that the theory of quantum gravity is a quantum theory of spacetime geometry with diffeo- 
morphism invariance (this legacy is discussed comprehensively in Rovelli's book [I122II ). To 
carry out the quantization procedure, one first casts general relativity into the Hamiltonian for- 
malism as a diff"eomorphism invariant Yang-Mills gauge field theory with a compact internal 
gauge group. Thus the construction of loop quantum gravity can also be applied to all back- 
ground independent gauge field theories. One can therefore claim that the theory can also be 
called as a background independent quantum gauge field theory. 

All classical fields theories, other than the gravitational field, are defined on a fixed space- 
time, which provides a foundation to the perturbative Fock space quantization. However general 
relativity is only defined on a manifold and hence is a background independent classical field 
theory, since gravity itself is the background. So the situation for gravity is much different 
from other fields by construction [ |122i |. namely gravity is not only the background stage, but 
also the dynamical actress. Such a double character for gravity leads to many difficulties in the 
understanding of general relativity and its quantization, since we cannot follow the strategy in 
ordinary quantum theory of matter fields. However, an amazing result in loop quantum gravity 
is that the background independent programme can even help us to avoid the difliculties in or- 
dinary quantum field theory. In perturbative quantum field theory in curved spacetime, the defi- 
nition of some basic physical quantities, such as the expectation value of energy-momentum, is 
ambiguous and it is difficult to calculate the back-reaction of quantum fields to the background 
spacetime [157]. One could speculate that the difliculty is related to the fact that the present 
formulation of quantum field theories is background dependent. For instance, the vacuum state 
of a quantum field is closely related to spacetime structure, which plays an essential role in the 
description of quantum field theory in curved spacetime and their renormalization procedures. 
However, if the quantization programme is by construction background independent and non- 
perturbative, it may be possible to solve the problems fundamentally. In loop quantum gravity 
there is no assumption of a priori background metric at all and the gravitational field and matter 
fields are coupled and fluctuating naturally with respect to each other on a common manifold. 

In the following sections, we will review pedagogically the basic construction of a com- 
pletely new, background independent quantum field theory, which is completely different from 
the known quantum field theory. For completeness and accuracy, we will use detailed math- 
ematical terminology. However, for simplicity, we will skip the complicated proofs of many 
important statements. One may find the missing details in the references cited. Thus our review 
will not be comprehensive. We refer to Ref. H 14611 and ||90l for a more detailed exploration, 
Refs. IfZDI and [I148II for more advanced topics. It turns out that in the framework of loop 
quantum gravity all theoretical inconsistencies introduced in the previous section are likely to 
be cured. More precisely, one will see that there is no UV divergence in quantum fields of 
matter if they are coupled with gravity in the background independent approach. Also recent 
works show that the singularities in general relativity can be smeared out in symmetry-reduced 
models [|43llll01||ll48l . The crucial point is that gravity and matter are coupled and consistently 
quantized non-perturbatively so that the problems of classical gravity and quantum matter in- 
consistency disappear. 
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2 Classical Framework 



2.1 Lagrangian Formalism 

In order to canonically quantize classical gravity, a Hamiltonian analysis has to be performed 
to obtain a canonical formalism of the classical theory suitable to be represented on a Hilbert 
space. A well known canonical formalism of general relativity is the ADM formalism (ge- 
ometrodynamics) derived from the Einstein-Hilbert actionl T36l[|971 . which has been problem- 
atic to cast into a quantum theory rigorously. Another well-known action of general relativity 
is the Palatini formalism, where the tetrad and the connection are regarded as independent dy- 
namical variables. However, the dynamics of the Palatini action has similar difficulties at the 
time of quantization as the dynamics derived from the Einstein-Hilbert action [ffl llSTl . In 1986, 
Ashtekar presented a formalism of true connection dynamics for general relativity with a rela- 
tively simple Hamiltonian constraint, and thus opened the door to apply quantization techniques 
from gauge field theory Il2lll3l lll231l . However, a drawback of that formalism is that the canon- 
ical variables are complex, which need the implementation of complicated reality conditions 
if one is to represent real general relativity. Moreover, the quantization based on the complex 
connection could not be carried out rigorously, since the internal gauge group is noncompact. 
In 1995, Barbero modified the Ashtekar new variables to give a system of real canonical vari- 
ables for dynamical theory of connections [34]. Then Hoist constructed a generalized Palatini 
action to support Barbero's real connection dynamics [|93l . Although there is a free parameter 
(Barbero-Immirzi parameter yS) in the generalized Palatini action and the Hamiltonian constraint 
is more complicated than the Ashtekar one, the generalized Palatini Hamiltonian with the real 
connections is widely used by loop theorists for the quantization programme. All the following 
analysis is based on the generalized Palatini formalism. 

Consider a 4-manifold, M, on which the basic dynamical variables in the generalized Pala- 
tini framework are a tetrad e" and an so{l,3)-va\ued connection oj^'''(not necessarily torsion- 
free), where the capital Latin indices /, 7, ... refer to the internal 5 0(1, 3) group and the Greek 
indices a, /3, ... denote spacetime indices. A tensor with both spacetime indices and internal 
indices is named as a generalized tensor. The internal space is equipped with a Minkowskian 
metric rjjj (of signature -, -I-) fixed once for all, such that the spacetime metric reads: 

The generalized Palatini action in which we are interested is given by [|201l : 

5p[4, = 1 £ d'x{e)e",(fij{aj/ + ^6"^,Q„/^), (2) 

where e is the square root of the determinant of the metric ga^, is the internal Levi-Civita 
symbol, ;S is the real Barbero-Immirzi parameter, and the 5o(l, 3)-valued curvature 2-form Q.^^^ 
of the connection co'/ reads: 

here Da denote the so{\,3) generalized covariant derivative with respect to oj^l^ acting on both 
spacetime and internal indices. Note that the generalized Palatini action returns to the Palatini 
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action when ^ = and gives the (anti)self-dual Ashtekar formalism when one sets ^ = +i. 
Moreover, besides spacetime diffeomorphism transformations, the action is also invariant under 
internal 5 0(1, 3) rotations: 

(e, oj) ^ (e', oj') = (b'^e, b-^cob + b'^db), 

for any 50(1,3) valued function b on M. The gravitational field equations are obtained by 
varying this action with respect to e" and oj^^ . We first study the variation with respect to the 
connection oj^^ One has 

by the definition of co variant generalized derivative Da- Note that dco^^ is a Lorentz covariant 
generalized tensor field since it is the difference between two Lorentz connections [|107lini04|| . 
One thus obtains 

5Sp = ^^£d\(e)e'^4(6nj' + ^e"^,6n,^^') 

where we have used the fact that Da A" = da A" for all vector density A" of weight +1 and 
neglected the surface term. Then it gives the equation of motion: 

where Ty"^^'' is the spacetime Levi-Civita symbol. This equation leads to the torsion-free Cartan's 
first equation: 

which means that the connection toj'' is the unique torsion-free Levi-Civita spin connection 
compatible with the tetrad e". As a result, the second term in the action Q can be calculated 
as: 

(e)e'^/je"'''^Cia^KL = n'"'''R./iy6, 

which is exactly vanishing, because of the symmetric properties of Riemann tensor. So the 
generalized Palatini action returns to the Palatini action, which will certainly give the Einstein 
field equation. 

2.2 Hamiltonian Formalism 

To carry out the Hamiltonian analysis of action Q, suppose the spacetime M is topologically 
2 X R for some 3-dimensional compact manifold E without boundary. We introduce a foliation 
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parameterized by a smooth function t and a time-evolution vector field f such that f{dt)a = 1 
in M, where t" can be decomposed with respect to the unit normal vector n" of S as: 

f = Nn" + N", (3) 

here A'^ is called the lapse function and A'^'^ the shift vector [|156ll[|97ll . The internal normal vector 
is defined as nj = n^e". It is convenient to carry out a partial gauge fixing, i.e., fix a internal 
constant vector field n' with rjjjn'n-' = -1. Note that the gauge fixing puts no restriction on 
the real dynamic^. Then the internal vector space V is 3+1 decomposed with a 3-dimensional 
subspace W orthogonal to n', which will be the internal space on S. With respect to the internal 
normal and spacetime normal n", the internal and spacetime projection maps are denoted by 
q' and q'^ respectively, where we use i, j, k, ...to denote the 3-dimensional internal space indices 
and a, b, c, ... to denote the indices of space 2. Then an internal reduced metric djj and a reduced 
spatial metric on Z, qab, are obtained by these two projection maps. The two metrics are related 
by: 

qab = Sije'^el, (4) 

where the orthonormal co-triad on S is defined by := e'^q',q". Now the internal gauge group 
50(1,3) is reduced to its subgroup S0{3) which leaves invariant. Finally, two Levi-Civita 
symbols are obtained respectively as 

I J K L 
^ijk '■= liljlk'^ ^LIJK, 

n := qlcf^qyri 

—abc " —iioipy 

where the internal Levi-Civita symbol 6,;^ is an isomorphism of Lie algebra 50(3). Using the 
connection 1-form (x)l\ one can defined two 5o(3)-valued 1-form on E: 

K ■■= q\q>"nj, 

where F is a spin connection on S and K will be related to the extrinsic curvature of S on shell. 
After the 3-1-1 decomposition and the Legendre transformation, action Q can be expressed as 



Sp= [dt [ d'x[P'J£A'a-not(Ai„P'j,A',N,N')], (5) 
Jr Js 

from which the symplectic structure on the classical phase space is obtained as 

{jV^X^XP'jiy)] := 6)6l6\x,y), (6) 
where the configuration and conjugate momentum are defined respectively by: 

Ai := 



'However, there are some arguments that such a gauge fixing is a non-natural way to break the internal Lorentz 
symmetry (see e.g. liOU ). 
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here det^ is the determinant of the 3-metric qab on 2 and hence det^ = (/c/3)^ det P. In the 
definition of the configuration variable A^, we should emphasize that Y'^ is restricted to be the 
unique torsion free 5'o(3)-valued spin connection compatible with the triad e". This conclusion 
is obtained by solving a second class constraint in the Hamiltonian analysis [|93l . In the Hamil- 
tonian formalism, one starts with the fields (A^,, P^). Then neither the basic dynamical variables 
nor their Poisson brackets depend on the Barbero-Immirzi parameter [5. Hence, for the case 
of pure gravitational field, the dynamical theories with diff"erent jS are related by a canonical 
transformation. However, as we will see, the spectrum of geometric operators are modified by 
different value of /3, and the non-perturbative calculation of black hole entropy is compatible 
with Bekenstein-Hawking's formula only for a specific value of [68J. In addition, it is argued 
that the Barbero-Immerzi parameter /3 may lead to observable eff"ects in principle when the grav- 
itational field is coupled with Fermions [|112| . In the decomposed action ([5]), the Hamiltonian 
density 'Htot is a linear combination of constraints: 



not = A'Gi + N"Ca + NC, 
N an 

Hamiltonian are expressed as GUI : 

■= BP" 
1+0^ 



where A' = -^e'jj^coj'', N" and are Lagrange multipliers. The three kinds of constraints in the 
^pressed as GUI : 



Gi -- 








c -- 


K0- 



2 VI det^l 



pa 

+ K{\ + 0')da{^^)G\ (7) 

V det^l 

where the configuration variable AJ, performs as a 5o(3)- valued connection on S and F^^^ is the 
5o(3)-valued curvature 2-form of with the well-known expression: 

:= 2D,X^ = daK - dtK + e'jAiK- 

In any dynamical system with constraints, the constraint analysis is essentially important be- 
cause they reflect the gauge invariance of the system. From the above three kinds of constraints 
of general relativity, one can know the gauge invariance of the theory. The Gauss constraint 
Gi = has crucial importance in formulating the general relativity into a dynamical theory of 
connections. The corresponding smeared constraint function, ^(A) := J^d^xAXx)Gi(x), gener- 
ates a transformation on the phase space as: 

{PlixXgiA)} = ei/A^(x)Pt(x), 

which are just the infinitesimal versions of the following gauge transformation for the so(3)- 
valued connection 1-form A and internal rotation for the 5o(3)-valued densitized vector field P 
respectively: 

(A„ P') ^ ig-'A„g + g-\dg),, g~'F'g). 
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To display the meaning of the vector constraint Ca = 0, one may consider the smeared constraint 
function: 



It generates the infinitesimal spatial diffeomorphism by the vector field A'^'' on E as: 

The smeared scalar constraint is weakly equivalent to the following function, which is re- 
expressed for quantization purpose as 



I' 



.S(A^) := I d'xN(x)C(x) 

papb 

i 

Vldet^l 



= '± fdhN-^ie^Ft,- 2(1+ fi'WlKi,]- (8) 
^ Js v deta 



It generates the infinitesimal time evolution off S. The constraints algebra, i.e., the Poisson 
brackets between these constraints, play a crucial role in the quantization programme. It can be 
shown that the constraints algebra of ^ has the following form: 

{^(A), ^(A')} = ^([A, A']), 

{^(A), ^(A^)} = -gu^A), 

mA),sm = 0, 

mA^), ^(A^')} = ^([A^, A^']), 
{^(A^), S(M)} = -SU^M), 
{SiN), SiM)} = -"ViiNdhM - MdbN)q"^) 

-g{{Nd,M -Mdt,N)q"'Aa)) 

I det q\ 

where |det^|^"'' = k^0^P"P^.6'K Hence the constraints algebra is closed under the Poisson 
brackets, i.e., the constraints are all first class. Note that the evolution of constraints is consistent 
since the Hamiltonian H = fj. d^x^tot is a linear combination of the constraints functions. The 
evolution equations of the basic canonical pair read 

= [K, H], Zri = [PI H}. 

Together with the three constraint equations, they are completely equivalent to the Einstein 
field equations. Thus general relativity is cast as a dynamical theory of connections with a 
compact structure group. Before finishing the discussion of this section, several issues should 
be emphasized. 
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• Canonical Transformation Viewpoint 

The above construction can be reformulated in the language of canonical transformations, 
since the phase space of connection dynamics is the same as that of triad geometrodynam- 
ics. In the triad formalism the basic conjugate pair consists of densitized triad E" = /SP'^ 
and "extrinsic curvature" K'^. The Hamiltonian and constraints read 

Hot = A'G'i+N"Ca+NC 
G; = e^'KiEt, (10) 

Ca = E'^^[aKiy (11) 



1 1 



C = [-1 det + E\"EfKiKl], (12) 



Vldet^l 
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where is the S 0(3) generalized derivative operator compatible with triad e° and R is 
the scalar curvature with respect to it. Since E" is a vector density of weight one, we have 

V^Et = daE^ + e,/riEt = 0. 
One can construct the desired Gauss law by 

G,- := -^V.Ef + G;, 

which is weakly zero by construction. This motivates us to define the connection A" = 
Y[^+/3K'^j. Moreover, the transformation from the pair (E", Kj^) to (P^ A^) can be proved to 
be a canonical transformation [13?ll[|146ll , i.e., the Poisson algebra of the basic dynamical 
variables is preserved under the transformation: 

E'l P^= E^/J3 

4 ^ Ai = ri+/3Ki, 



{Ff(x), Aiiy)} = {E'^{x),Ki(y)} = 5l5\6{x,y), 
{A[{x\Aiiy)} = {K'^ixXKiiy)} = 0, 
{P^'(x), P){y)] = {E';(x\E'(y)} = 0. 



• The Preparation for Quantization 

The advantage of a dynamical theory of connections is that it is convenient to be quan- 
tized in a background independent fashion. In the following procedure of quantization, 
the quantum algebra of the elementary observables will be generated by holonomy, i.e., 
connection smeared on a curve, and electric flux, i.e., a densitized triad smeared on a 2- 
surface. So no information about a background geometry is needed to define the quantum 
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algebra. In the remainder of the thesis, in order to incorporate also spinors, we will en- 
large the internal gauge group to be SU{2). This does not damage the prior constructions 
because the Lie algebra of S U(2) is the same as that of 50(3). Due to the well-known 
nice properties of compact Lie group S U{2), such as the Haar measure and Peter- Weyl 
theorem, one can obtain the background independent representation of the quantum alge- 
bra and the spin-network decomposition of the kinematic Hilbert space. 

• Analysis on Constraint Algebra 

The classical constraint algebra ^ is an infinite dimensional Poisson algebra. Unfortu- 
nately, it is not a Lie algebra since the Poisson bracket between two scalar constraints has 
structure function depending on dynamical variables. This causes problems when solving 
the constraints quantum mechanically. On the other hand, one can see from Eq.® that 
the algebra generated by Gauss constraints forms not only a subalgebra but also a 2-side 
ideal in the full constraint algebra. Thus one can first solve the Gauss constraints indepen- 
dently. It is convenient to find the quotient algebra with respect to the Gauss constraint 
subalgebra as 



which plays a crucial role in solving the constraints quantum mechanically. But the subal- 
gebra generated by the difFeomorphism constraints can not form an ideal. Hence the pro- 
cedures of solving the diffeomorphism constraints and solving Hamiltonian constraints 
are entangled with each other. This leads to certain ambiguity in the construction of a 
Hamiltonian constraint operator P 134||||149 i Fortunately, the Master Constraint Project 
addresses the above two problems as a whole by introducing a new classical constraint al- 
gebra II149II . The new algebra is a Lie algebra where the difFeomorphism constraints form 
a 2-side ideal. We will come back to this point in the discussion on quantum dynamics of 
loop quantum gravity. 



{^(A^), SiM)] 
{S(N), S(M)} 
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3 Foundations of Loop Quantum Gravity 



In this chapter, we will begin to quantize the above classical dynamics of connections as a 
background independent quantum field theory. The main purpose of the chapter is to construct 
a suitable kinematical Hilbert space 'Hkm for the representation of quantum observables. In the 
following discussion, we formulate the construction in the language of algebraic quantum field 
theory [|85l . It should be emphasized that the following constructions can be generalized to all 
background independent non-perturbative gauge field theories with compact gauge groups. 

3.1 General Programme for Algebraic Quantization 

In the strategy of loop quantum gravity, a canonical programme is performed to quantize gen- 
eral relativity, which has been cast into a diffeomorphism invariant gauge field theory, or more 
generally, a dynamical system with constraints. The following is a summary for a general pro- 
cedure to quantize a dynamical system with first class constraints [il33il [i4il. 

• Algebra of Classical Elementary Observables 

One starts with the classical phase space {A\, {, }) and R (R can be countable infinit>@ ) 
first-class constraints C,.(r = 1...R) such that {€,-, Cv) = l^f^^f^/Ct, where is generally 
a function on phase space, namely, the structure function of Poisson algebra. The algebra 
of classical elementary observables ^ is defined as: 

Definition 3.1.1: The algebra of classical elementary observables ^ is a collection of 
functions f{m), m e Ai on the phase space such that 

(1) /(m) 6 should separate the points of Ai, i.e., for any m m', there exists f{m) 6 ^, 
such that fim) f(m'); {analogy to the p and q in Ai = T*R.) 

(2) /(m), f'On) 6 => {f(m), f'(m)} e ^ (closed under Poisson bracket); 

(3) /(m) e => /(m) 6 'ip (closed under complex conjugation). 

So forms a sub *-Poisson algebra of C'^(A\). In the case of AI = T*R, is gener- 
ated by the conjugate pair (q, p) with {q, p] = \. 

• Quantum Algebra of Elementary Observables 

Given the algebra of classical elementary observables S^, a quantum algebra of elementary 
observables can be constructed as follows. Consider the formal finite sequences of classi- 
cal observable (/i.../„) with 6 ^. Then the operations of multiplication and involution 
are defined as 

(/l > •••5 ■ (/l J •••5 /m) •~ ifl: •■■^ fm fl^ •■■^ fin)'' 
(fl,--,fn)* '■= (fn,--;fl)- 

One can define the formal sum of different sequences with different number of elements. 
Then the general element of the newly constructed free * -algebra FQl^i) of is formally 

^This includes the case of field theory with infinite many degree of freedom, since one can introduce the 
expression C,,,;^ - cPx^n{x)C^{x), where {0n(x))JJ!^j forms a system of basis in L^(E, cPx). 
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expressed as Yjk=i(f\'^^ where /^'^ 6 Consider the elements of the form (se- 

quences consisting of only one element) 

(/ + /')-(/)-(/'), izf)-zif), [(/),(/')]- W,/'}), 
where z e C is a complex number, and the canonical commutation bracket is defined as 

[(/),(/')] :=(/)•(/')-(/')■(/)• 

A 2-side ideal 3 of F^^) can be constructed from these element, and is preserved by the 
action of involution *. One thus obtains 

Definition 3.1.2: The quantum algebra % of elementary observables is defined to be the 
quotient *-algebra FC^)/^. 

Note that the motivation to construct a quantum algebra of elementary observables is 
to avoid the problem of operators ordering in quantum theory so that the quantum algebra 
% can be represented on a Hilbert space without ordering ambiguities. 

• Representation of Quantum Algebra 

In order to obtain a quantum theory, we need to quantize the classical observables in the 
dynamical system. The, so called, quantization is nothing but a * -representation majj^ n 
from the quantum algebra of elementary observable 91 to the collection of linear opera- 
tors JICH) on a Hilbert Space 'H. At the level of quantum mechanics, the well-known 
Stone- Von Neumann Theorem concludes that in quantum mechanics, there is only one 
strongly continuous, irreducible, unitary representation of the Weyl algebra, up to uni- 
tary equivalence (see, for example, Ref. FllBH ). However, the conclusion of Stone- Von 
Neumann cannot be generalized to the quantum field theory because the latter has infinite 
many degrees of freedom (for detail, see, for example [157J). In quantum field theory, 
a representation can be constructed by GNS(Gerfand-Naimark-Segal)-construction for a 
quantum algebra of elementary observable 91, which is a unital *-algebra actually. The 
GNS-construction for the representation of quantum algebra 91 is briefly summarized as 
follows. 

Definition 3.1.3: Given a positive linear functional (a state) co on % the null space 
JfJ^ 6 91 with respect to co is defined a* ^R^ := {a e 9I|a>(a* • a) = 0}, which is a left ideal 
in 91. Then a quotient map can be defined as [.]: 91 ^/yioj,' a i-^ [a] := {a + b\b e ytai}- 

map n: % ^ -CCH) is a * -representation if and only if (1) there exists a dense subspace D of 1-1 contained 
in f^aei\.{D{n(a)) n D{ji(a*))\ where D{ji{a)) is the domain of the operator ji{a) and (2) for every a,b e'H and A e C 
the following conditions are satisfied in £), 

7T(a + b) - 7r(a) + n(b), n(Aa) - An(a), 
n(a ■ b) = n(a)n{b), nia*) - n{a)\ 

Note that JLCH) fails to be an algebra because the domains of unbounded operators cannot be the whole Hilbert 
space. However, the collection of bounded operators on any Hilbert space is really a *-algebra. 
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The GNS-representation for '31 with respect to co is a ^-representation map: nc^: '31 
X('K^), where 'Ha, '■= C^/^cj) iind {.) denotes the completion with respect to the naturally 
equipped well-defined inner product 

< MIW >'H^'-= <^(«* • b) 

on "K^. This representation map is defined by 

n^{a){b] := [a • Zj], V a 6 « and M 6 "K^, 

where n^jia) is an unbounded operator in general. Moreover, GNS-representation is a 
cyclic representation, i.e., 3 Q-^j e 'Hoj, such that {{n{a)Q.J\a 6 %}) = IHoj and Q-^j is called 
a cyclic vector in the representation space. In fact Q.^) := [1] is a cyclic vector in 'Ha, and 
{[na,{a)Q.J\a e ^}) = "K^. As a result, the positive linear functional with which we begin 
can be expressed as 

a)(a) =< Q.j7Taj(.a)Q.aj >'K^ • 

Thus a positive linear functional on % is equivalent to a cyclic representation of%, which 
is a triple ("K^jTr^, ^a,)- Moreover, every non-degenerate representation is an orthogonal 
direct sum of cyclic representations {for proof, see ||5SI ) . 

So the kinematical Hilbert space 'Hkin = 'Ha, for the system with constrains can be ob- 
tained by GNS-construction. In the case that there are gauge symmetries in our dy- 
namical system, supposing that there is a gauge group G acting on % by automorphisms 
ag : % ^ %, V ^ 6 G, it is preferred to construct a gauge invariant representation of '31. 
So we require the positive linear functional oj on ?I to be gauge invariant, i.e., a)oag = oj. 
Then the group G is represented on the Hilbert space 'Ha, as: 

U{g)na,{a)Q.aj = na,{ag{a))Qaj, 

and such a representation is a unitary representation of G. In loop quantum gravity, it is 
crucial to construct an internal gauge invariant and diffeomorphism invariant representa- 
tion for the quantum algebra of elementary observables. 

• Implementation of the Constraints 

In the Dirac quantization programme for a system with constraints, the constraints should 
be quantized as some operators in a kinematical Hilbert space 'Hkin- One then solves 
them at quantum level to get a physical Hilbert space 'Hphys, that is, to find a quantum 
analogy C, of the classical constraint formula Cr and to solve the general solution of 
the equation = 0. However, there are several problems in the construction of the 
constraint operator C, . 

(i) C, is in general not in^, so there is a factor ordering ambiguity in quantizing Cr to 
be an operator Cr. 

(ii) In quantum field theory, there are ultraviolet(UV) divergence problems in construct- 
ing operators. However, the UV divergence can be avoided in the background inde- 
pendent approach. 
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(iii) Sometimes, quantum anomalies may appear when there are structure functions in 
the Poisson algebra. Although classically we have {C^, C,} = 2f^j/„'Q, r,s,t = 
1, ...,R, where f^J is a function on phase space, quantum mechanically it is possible 
that [C,., Cs\ 4^ ifiLf^-^f^jC, due to the ordering ambiguity between f^J and C,. If one 
sets [C„ C,s\ = '^llf^^iijCt + CtfJ), for ^ satisfying C.^P = 0, we have 

[c„ c,]Y = - 2] = - 2][a, /,./]^. (13) 

(=1 (=1 

However, [Cf,/^/]*P are not necessary to equal to zero for all r,s,t = I...R. If 
not, the problem of quantum anomaly appears and the new quantum constraints 
[C,,/^/]*P = have to be imposed on physical quantum states, since the classical 
Poisson brackets {C^, C^} are weakly equal to zero on the constraint surface M c M. 
Thus too many constraints are imposed and the physical Hilbert space l-fphys would 
be too small. Therefore this is not a satisfactory solution and one needs to find a 
way to avoid the quantum anomalies. 

• Solving the Constraints and Physical Hilbert Space 

In general the original Dirac quantization approach can not be carried out directly, since 
there is usually no nontrivial 6 "K^.,,, such that C^^P = 0. This happens when the con- 
straint operator C,- has "generalized eigenfunctions" rather than eigenfunctions. One then 
develops the so-called Refined Algebraic Quantization Programme, where the solutions 
of the quantum constraint can be found in the algebraic dual space of a dense subset in 
'Hkin (see e.g. [l84l ). The quantum diffeomorphism constraint in loop quantum gravity is 
solved in this way. Another interesting way to solve the quantum constraints is the Mas- 
ter Constraint Approach proposed by Thiemann recently ni49L which seems especially 
suited to deal with the particular feature of the constraint algebra of general relativity. A 
master constraint is defined as M := ^Z^^.^jK^sC^Cr for some real positive matrix K^^. 
Classically one has M = if and only if = for all r = I...R. So quantum mechan- 
ically one may consider solving the Master Equation: M¥ = to obtain the physical 
Hilbert space l-fphys instead of solving C^*? = 0, V r = 1...R. Because the master con- 
straint M is classically positive, one has opportunities to implement it as a self-adjoint 
operator on T-lkin- If it is indeed the case and 'Hkin is separable, one can use the direct in- 
tegral representation of 9ikin associated with the self-adjoint operator M to obtain ^phys- 

'Hkin - J^diu(A)'Hl 
<0|¥>fa„ = £diu(A)<(!>\^>^e, (14) 

where is a so-called spectral measure and 'Hf is the (generalized) eigenspace of M with 
the eigenvalue A. The physical Hilbert space is then formally obtained as 'Hphys = *T^x=o 
with the induced physical inner product < | >'Kf_g 0- Now the issue of quantum anomaly 

One need to be careful for such a formal prescription, see the later discussion of master constraint or li63l . 
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is represented in terms of the size of l-iphys and the existence of sufficient numbers of 
semi-classical states. 

• Physical Observables 

We denote A\ as the original unconstrained phase space, M as the constraint surface, i.e., 
M := {m e M\Cr{rn) = 0, V r = I...R}, and M as the reduced phase space, i.e. the space 
of orbits for gauge transformations generated by all C^. The concept of Dirac observable 
is defined as the follows. 



Definition 3.1.4: 

( 1 ) A function O on M is called a weak Dirac observable if and only if the function de- 
pends only on points of M, i.e., {O, C^}|y^ = 0/or all r = I...R. For the quantum version, 
a self-adjoint operator O is a weak Dirac observable if and only if the operator can be 
well defined on the physical Hilbert space. 

(2) A function O on M is called a strong Dirac observable if and only if {0,Cr}\M = 
for all r = 1...R. For the quantum version, a self-adjoint operator O is a strong Dirac 
observable if and only if the operator can be defined on the kinematic Hilbert space "Hun 
and \p, Cr\ = in l-ikmfor all r = I...R. 

A physical observable is at least a weak Dirac observable. While Dirac observables have 
been found explicitly in symmetry reduced models, some even with an infinite number 
of degrees of freedom, it seems extremely difficult to find explicit expressions for them 
in full general relativity. Moreover the Hamiltonian is a linear combination of first-class 
constraints. So there is no dynamics in the reduced phase space, and the meaning of 
time evolution of the Dirac observables becomes subtle. However, using the concepts of 
partial and complete observables II121IIII115IIII122L a systematic method to get Dirac ob- 
servables can be developed, and the problem of time in such system with a Hamiltonian 
H = ll^^^PrCr may also be solved. 

Classically, let /(m) and {Tr{m)]^^^ be gauge non-invariant functions (partial observables) 
on phase space A1, such that Asr = {Cv, Tr} is a non-degenerate matrix on M. A system 
of classical weak Dirac observables (complete observables) FJ^, labelled by a collection 
of real parameters r = {r^jf^j can be constructed as 



F^ - V (^i--^ir---(T^-^i^)"'- y„. ^ 



where X,(/) := {Sf^jA-'C,, /} = {C„ /} . It can be verified that [X„ = and 

{Fj.^, CA\m = 0' for all r = I...R (for details see MB and [|62||). 

The partial observables {Tr(m)}f^^ may be regarded as clock variables, and is the time 
parameter for T^. The gauge is fixed by giving a system of functions {Tr(m)}f^^ and 
corresponding parameters {r^jf^j, namely, a section in M is selected by T^im) = t,. for 
each r, and Fj.j. is the value of / on the section. To solve the problem of dynamics, one 
assumes another set of canonical coordinates (Pi,-- -, Pn-r, Hi, • • •, H^; Qi,---, Qn-r, Ti,-- 
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•, Tr) by canonical transformations in the phase space (A1, { , }), where and are 
conjugate to 2^ and respectively. After solving the complete system of constraints 
{C,.(Pi, Qj, n,, Tt) = 0}^^j, the Hamiltonian H,. with respect to the time is obtained as 
Hr := TlriPi, Qj, Tt). Given a system of constants {(To)r}f=i, for an observable /(P,, Qj) 
depending only on and Qj, the physical dynamics is given by liMI [115011 : 

d 

where j is the physical Hamiltonian function generating the evolution with respect to 
Tr- Thus one has addressed the problem of time and dynamics as a result. 

• Semi-classical Analysis 

An important issue in the quantization is to check whether the quantum constraint oper- 
ators have correct classical limits. This has to be done by using the kinematical semi- 
classical states in •K|,„. Moreover, the physical Hilbert space l-iphys must contain enough 
semi-classical states to guarantee that the quantum theory one obtains can return to the 
classical theory when % ^ 0. The semi-classical states in a Hilbert space 'H should have 
the following properties. 

Definition 3.1.5: Given a class of observables S which comprises a subalgebra in the 
space X('K) of linear operators on the Hilbert space, a family of (pure) states {cOm}meM 
are said to be semi-classical with respect to S if and only if: 

(1) The observables in S should have correct semi-classical limit on semi-classical states 
and the fluctuations should be small, i.e., 

w,n{a)-aim) 
lim| 

n^o a{m) 

ojmia^) - <x),„iaf 
lim I 

for all a e S. 

(2) The set of cyclic vectors Q.^ related to co^ via the GNS -representation (n 1^,9^01,^0)) 
is dense in 9i. 




Seeking for semiclassical states are one of open issues of current research in loop quantum 
gravity. Recent original works focus on the construction of coherent states of loop quan- 
tum gravity in analogy with the coherent states for harmonic oscillator system II142IIII143II 

The above is the general programme to quantize a system with constraints. In the following sub- 
section, we will apply the programme to the theory of general relativity and restrict our view to 
the representation with the properties of background independence and spatial diffeomorphism 
invariance. 
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3.2 Quantum Configuration Space 

In quantum mechanics, the kinematical Hilbert space is L^(R^, d^x), where the simple is the 
classical configuration space of free particle which has finite degrees of freedom, and d^x is the 
Lebesgue measure on R^. In quantum field theory, it is expected that the kinematical Hilbert 
space is also the space on the configuration space of the field, which is infinite dimensional, 
with respect to some Borel measure naturally defined. However, it is often hard to define a 
concrete Borel measure on the classical configuration space, since the integral theory on infinite 
dimensional space is involved [57|. Thus the intuitive expectation should be modified, and the 
concept of quantum configuration space should be introduced as a suitable enlargement of the 
classical configuration space so that an infinite dimensional measure, often called cylindrical 
measure, can be well defined on it. The example of a scalar field can be found in the refer- 
ences f20]l E4l . For quantum gravity, it should be emphasized that the construction for quantum 
configuration space must be background independent. Fortunately, general relativity has been 
reformulated as a dynamical theory of SU (2) connections, which would be great helpful for our 
further development. 

The classical configuration space for gravitational field, which is denoted by Jl, is a col- 
lection of the 5M(2)-valued connection 1-form field smoothly distributed on E. The idea of the 
construction for quantum configuration is due to the concept of holonomy. 

Definition 3.2.1: Given a smooth S U(2) connection field A'^ and an analytic curve c with the 
parameter t e [0, 1] supported on a compact subset (compact support ) ofZ, the corresponding 
holonomy is defined by the solution of the parallel transport equation !\104\1 

jA{c,t) = -{A\eTi\A{c,t\ (15) 

with the initial value A{c, 0) = 1, where c" is the tangent vector of the curve and Tj e su(2) 
constitute an orthonormal basis with respect to the Killing-Cartan metric t](^,0 •= ~2Tr(^^), 
which satisfy [r,-, Tj] = ^'^ijTk and are fixed once for all. Thus the holonomy is an element in 
S U (2), which can be expressed as 

A(c)=!Pexp(- r [Aic'^n] dt), (16) 
Jo 

where A{c) e S U(2) and V is a path-ordering operator along the curve c (see the footnote at 
p382 in n04\l ). 

The definition can be well extended to the case of piecewise analytic curves via the relation: 

A(cioc2)=A(ci)A(c2), (17) 

where o stands for the composition of two curves. It is easy to see that a holonomy is invariant 
under the re-parametrization and is covariant under changing the orientation, i.e., 

A{c-')=A{crK (18) 
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So one can formulate the properties of holonomy in terms of the concept of the equivalent 
classes of curves. 

Definition 3.2.2: Two analytic curves c and c' are said to be equivalent if and only if they 
have the same source sic) {beginning point ) and the same target tic) {end point ), and the 
holonomies of the two curves are equal to each other, i.e., A(c) = A(c') VA 6 JH. A equivalent 
class of analytic curves is defined to be an edge, and a piecewise analytic path is an composition 
of edges. 

To summarize, the holonomy is actually defined on the set f of piecewise analytic paths with 
compact supports. The two properties (fTTl) and (fTSl) mean that each connection in J?l is a ho- 
momorphism from which is so-called a groupoid by definition il55 l|, to our compact gauge 
group SUil). Note that the internal gauge transformation and spatial diffeomorphism act co- 
variantly on a holonomy as 

A{e) ^ gitie))~^Aie)gisie)) and A(e) ^ o e), (19) 

for any 5 f/(2)-valued function gix) on E and spatial diffeomorphism (p. All above discussion 
is for classical smooth connections in The quantum configuration space for loop quantum 
gravity can be constructed by extending the concept of holonomy, since its definition does not 
depend on an extra background. One thus obtains the quantum configuration space Jl of loop 
quantum gravity as the following. 

Definition 3.2.3: The quantum configuration space ^ is a collection of all quantum connections 
A, which are algebraic homomorphism maps without any continuity assumption from the col- 
lection of piecewise analytic paths with compact supports, P, on D to the gauge group SUil), 
^:=Hom(^,5[/(2)fl 

Thus for any A e ^ and edge e in P, 

Aiex o e2) = A(ei)A(e2) and Aie'^) = Aie)-\ 

The transformations of quantum connections under internal gauge transformations and diffeo- 
morphisms are defined by Eq.U^. 

The above discussion on the smooth connections shows that the classical configuration space 
^ can be understood as a subset in the quantum configuration space Moreover, the Giles 
theorem |[82l shows precisely that a smooth connection can be recovered from its holonomies 
by varying the length and location of the paths. 

On the other hand, it was shown in H155IIII146II that the quantum configuration space J[ can 
be constructed via a projective limit technique and admits a natural defined topology. To make 
the discussion precise, we begin with a few definitions: 

Definition 3.2.4: 

^It is easy to see that the definition of J[ does not depend on the choice of local section in S f/(2)-bundle, since 
the internal gauge transformations leave J[. invariant. 
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1. A finite set {e\, e^] of edges is said to be independent if the edges e,- can only intersect 
each other at their sources s(ei) or targets t{ei). 

2. A finite graph is a collection of a finite set {ei,...,eN} of independent edges and their 
vertices, i.e. their sources s{ei) and targets t(ei). We denote by E{y) and y(y) respectively 
as the sets of independent edges and vertices of a given finite graph y. And Ny is the 
number of elements in E{y). 

3. A subgroupoid a(y) c f can be generated from y by identifying V{y) as the set of ob- 
jects and all e e E{y) together with their inverses and finite compositions as the set of 
homomorphisms. This kind of subgoupoid in V is called tame subgroupoid. a{y) is inde- 
pendent of the orientation ofy, so the graph y can be recovered from tame subgroupoid 
a up to the orientations on the edges. We will also denote by Na the number of elements 
in E{y) where y is recovered by the tame subgroupoid a. 

4. X denotes the set of all tame subgroupoids in V. 

One can equip a partial order relation < on xlE defined hy a < a' if and only if c is a sub- 
groupoid in a' . Obviously, for any two tame subgroupoids a = a(y) and a' = a{y') in X, there 
exists a" = a{y") 6 X such that a, a' < a", where y" = y (J y' . Define Xa = Hom(a, S U(2)) 
as the set of all homomorphisms from the subgroupoid a(y) to the group SU{2). Note that an 
element Aa e Xa (a = a(y)) is completely determined by the S U (2) group elements A(e) where 
e e E(y), so that one has a bijection A : Xa(y) SU{2)'^^, which induces a topology on Xa{y) 
such that /I is a topological homomorphism. For any pair a < a', one can define a surjective 
projection map Pa'a from Xa' to Xg, by restricting the domain of the map Aa' from a' to the sub- 
groupoid a, and these projections satisfy the consistency condition P^'a ° Pa" a' = Pa" a- Thus 
a projective family {Xa, Pa'a]a<a' is obtained by above constructions. Then the projective limit 
limaiXa) is naturally obtained. 

Definition 3.2.5: The projective limit \ima(Xa) of the projective family {Xa,Pa'a}a<a' is a subset 
of the direct product space Xoo '■= Y\ae£ ^a defined by 

lim(Xa) := {{Aa}ae£\Pa'aAa' = Aa, a < Or'}. 

a 

Note that the projection Pa'a is surjective and continuous with respect to the topology of Xa. 
One can equip the direct product space Xoo with the so-called Tychonov topology. Since any Xa 
is a compact Hausdorff space, by Tychonov theorem Xoo is also a compact Hausdorff space. One 
then can prove that the projective limit, Yima{Xa), is a closed subset in Xoo and hence a compact 
Hausdorff space with respect to the topology induced from Xoo. At last, one can find the relation 
between the projective limit and the prior constructed quantum configuration space J{. As one 
might expect, there is a bijection O between J?l and limQ,(XQ,) [I146II : 

O: ^ lim(X^); 

a 

A 1-^ {A\a}ae£, 

partial order on X is a relation, which is reflective (a < a), symmetric {a < a' , a' < a ^ a' - a) and 
transitive (a < a', a' < a" ^ a' < a"). Note that not all pairs in X need to have a relation. 
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where A\a means the restriction of the domain of the map A e Jl = Hom(^, S U(2)). As a result, 
the quantum configuration space is identified with the projective Umit space and hence can be 
equipped with the topology. In conclusion, the quantum configuration space J{. is constructed 
to be a compact HausdorfF topological space. 

3.3 Cylindrical Functions on Quantum Configuration Space 

Given the projective family [Xa, Pa'a}a<a', the cylindrical function on its projective limit is 
well defined as follows. 

Definition 3.3.1: Let C{Xa) be the set of all continuous complex functions on Xq,, two func- 
tions fa G C{Xa) and fa' e C{Xa') are said to be equivalent or cylindrically consistent, denoted 
by fa ~ fa', if and only ifP*a»afa - P*a"a'fa'' ^ ^' ^' ' where P*„^ denotes the pullback map 
induced from Pa" a- Then the space Cyl(Ji) of cylindrical functions on the projective limits is 
defined to be the space of equivalent classes [/], i.e., 

Cym):=[Ua C(Xa)y~. 

One then can easily prove the following proposition by definition. 
Proposition 3.3.1: 

All continuous functions fa on Xa are automatically cylindrical since each of them can generate 
a equivalent class {fa\ via the pullback map P*a,af^^ ^' ^ ^' '^^'^ ^^e dependence of P*^,^fa 
on the groups associated to the edges in or' but not in a is trivial, i.e., by the definition of the 
pull back map, 

(Pl'afaXAiei), ...,A(e^J, ...,AieNj) = fa(A(e,), A(e^J), (20) 

where Na denotes the number of independent edges in the graph recovered from the groupoid a. 
On the other hand, by definition, given a cylindrical function f e Cyl{^) there exists a suitable 
groupoid a such that f = [fa], so one can identify f with fa. Moreover, given two cylindrical 
functions f, f e Cyl(M), by definition of cylindrical functions and the property of projection 
map, there exists a common groupoid a and fa, fa^ C(Xa) such that f = [fa] and f = [fa]. 

Let /, /' e Cyl(Jl), there exists graph a such that / = [fa], and /' = [f^], then the following 
operations are well defined 

/ + /':= [fa + fa], ff ■■= [fafa], Zf := [zfa], f := [/.], 

where z e C and / denotes complex conjugate. So we construct Cyl(Jl) as an Abelian *-algebra. 
In addition, there is a unital element in the algebra because Cyl(Ji) contains constant functions. 
Moreover, we can well define the sup-norm for / = [fa] by 

11/11 := sup \fa(Aa)\, (21) 
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which satisfies the C* property ||//|| = Then Cyl{^) is a unital Abelian C* -algebra, after 
the completion with respect to the norm. 

From the theory of C* -algebra, it is known that a unital Abelian C* -algebra is identical to the 
space of continuous functions on its spectrum space via an isometric isomorphism, the so-called 
Gel'fand transformation (see e.g. 1114611 '). So we have the following theorem [fTTllfTSl . which 
finishes this section. 

Theorem 3.3.1: 

(1) The space Cyl{Jl) has a structure of a unital Abelian C* -algebra after completion with re- 
spect to the sup-norm. 

(2) Quantum configuration space 3\ is the spectrum space of completed Cyl{fK) such that 
CylifK) is identical to the space C{fK) of continuous functions on 3\. 

3.4 Loop Quantum Kinematics 

In analogy with the quantization procedure of section 3.1, in this subsection we would like to 
perform the background-independent construction of algebraic quantum field theory for gen- 
eral relativity. First we construct the algebra of classical observables. Taking account of the 
future quantum analogs, we define the algebra of classical observables ^ as the Poission *- 
subalgebra generated by the functions of holonomies (cylindrical functions) and the fluxes of 
triad fields smeared on some 2-surface. Namely, one can define the classical algebra in analogy 
with geometric quantization in finite dimensional phase space case by the so-called classical 
Ashtekar-Corichi-Zapata holonomy-flux *-algebra as the following [1961. 

Definition 3.4.1 

The classical Ashtekar-Corichi-Zapata holonomy-flux ^-algebra is defined to be a vector space 
^Acz := Cyl(Ji) X "V^^Ji), where "V^iJl) is the vector space of algebraic vector fields spanned 
by the vector fields if/Yf{S) if/ e Cyl{M), and their commutators, here the smeared fiux vector 
field YfiS ) is defined by acting on any cylindrical function: 



for any su(2)-valued function f with compact supports on S and if/ are cylindrical functions on 
M. We equip ^acz with the structure of an *-Lie algebra by: 
(1) Lie bracket { , } : ^acz x '^acz ''^acz is defined by 



(2) Involution: p p p e ^acz is defined by complex conjugate of cylindrical functions 
and vector fields, i.e., p := (if/, Y)'i p = (if/, Y) £ ^acz, where Yif/ := Yif/. 




{(if,, Y), (if/', Y')} := (Y o if,' - r o ifr, [Y, Y']), 



for all (if/, Y), (if/', Y') e ^acz with if/, if/' e Cyl(Jl) and Y, Y' e ^'^(^). 
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(3) ^Acz admits a natural action ofCyl{J{) by 
which gives ^acz a module structure. 

Note that the action of flux vector field 7/(5 ) on can be expressed explicitly on any cylindrical 
function if/y e C^{Xa{y)) via a suitable regularization 1114611 : 



's ' 

...... s-"'' ^^(^)- 

Sne,s(e) 



eeE(y) ^ OA{e)„ 
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where A{e)mn is the 81/(2) matrix element of the holonomy along the edge e, X':'^'^'^ is the 
left(right) invariant vector field L^'^'\R^'^'^) of the group associated with the edge e if v is the 
source(target) of edge e by definition: 

L(^-V(A(e)) := ^l=oHA{e)cxp(tTd), 
dt 

i?(^'V(A(e)) := i,^o<A(exp(-?T,)A(e)), 
dt 

and 

{0, if e n 5 = 0, or e lies in S ; 
1, if e lies above S and e OS = p; 
-1, if e lies below S and e OS = p. 

Since the surface S is oriented with normal n^, "above" means n^e^lp > 0, and "below" means 
na^"\p < O7 where e"\p is the tangent vector of e at p. And one should consider e r\ S contained 
in the set V{y) and some edges are written as the union of elementary edges which either lie in 
S , or intersect S at their source or target. On the other hand, from the commutation relations 
for the left(right) invariant vector fields, one can see that the commutators between flux vector 
fields do not necessarily vanish when S Ci S' (/). This unusual property is the classical origin 
of the non-commutativity of quantum Riemannian structures [j23l. 

The classical Ashtekar-Corichi-Zapata holonomy-flux *-algebra serves as a classical alge- 
bra of elementary observables in our dynamical system of gauge fields. Then one can construct 
the quantum algebra of elementary observables from ^acz in analogy with Definition 3.1.2. 

Definition 3.4.2[|96l 
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The abstract free algebra F{^acz) of the classical *-algebra is defined by the formal direct sum 
of finite sequences of classical observables {pi, ...,p„) with pk e ^acz. where the operations of 
multiplication and involution are defined as 

ipi,...,Pn)-(p\,...,p'J := ipi,...,Pn,p'i,...,p'J, 

iPu-,PnT ■= (Pn,-,Pl)- 

A 2-sided ideal 3 can be generated by the following elements, 

(p + p') - (p) - (p), (zp) - zip), 
m,ip')]-ih({p,p'}), 
((<A,0),;?)-(iAop), 

where the canonical commutation bracket is defined by 

[ip),(p')]:=ip)-ip')-(p')-ip). 

Note that the ideal 3 is preserved by the involution *, and the last set of generators in the ideal 
3 cancels the overcompleteness generated from the module structure of^Acz /E!/- 
The quantum holonomy-fiux *-algebra is defined by the quotient *-algebra ?I = FC^acz)/^, 
which contains the unital element 1 := ((1,0)). Note that a sup-norm has been defined by 
Eq.i2T]i for the Abelian sub-*-algebra Cyl(Jl) in % 

For simplicity, we denote the one element sequences (equivalence classes) {(i/f, 0)) and ((0, Y)) 
'i iff e Cyl{Jl), Y 6 "V^^iJl) in '21 by iff and Y respectively, where the "hat" denotes the equiv- 
alence class with respect to the quotient. In particular, for all cylindrical functions and flux 
vector fields Yf(S ), 

{fr* = ^ and YfiSy = 7/5). 

It can be seen that the free algebra F(^acz) is simplified a great deal after the quotient, and every 
element of the quantum algebra % can be written as a finite linear combination of elements of 
the form 

-A, 

• Yf^,{S yti) • Yf^^(S k2) ■ ■■■ • Yff^^iS kk). 

Moreover, given a cylindrical function iff and a flux vector field Yf(S), one has the relation from 
the commutation relation: 

7/(5 ) • (A = ihYj{S)if/ + (/r • Yf{S ) . (22) 
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Then the kinematical Hilbert space 'Hun can be obtained properly via the GNS -construction for 
unital *-algebra '21 in the same way as in Definition 3.1.3. By the GNS-construction, a positive 
linear functional, i.e. a state ojkim on % defines a cyclic representation ("Kfavj, nkin, ^kin) for '21. In 
our case of quantum holonomy-flux *-algebra, the state with both Yang-Mills gauge invariance 
and diffeomorphism invariance is defined for any ij/y e Cyl{^) and non- vanishing flux vector 
field 7/(5) G^^C;^) as [|96l: 

(Okini^fy) ■= I F[ dlXH(A(e))lfry({A(e)}eeE(y)), 

JsUilfr ,^E^y^ 

aJkin(a-Yf{S)):=0, \/ae% 

where dfin is the Haar measure on the compact group SU(2) and Ny is the number of elements 
in E{y). This ojkm is called Ashtekar-Isham-Lewandowski state. The null space yikin £ ^ with 
respect to totm is defined as 3\i„ := {a e %\a)kinia* ■ a) = 0}, which is a left ideal. Then a quotient 
map can be defined as: 

[.] : « ^ 

a 1-^ [a] := {a + b\b e "^km}- 

The GNS-representation for '21 with respect to tOkm is a representation map: TT/t,,,: '21 -Li^kin) 
such that nkinia ■ b) = nkm{a)nkin(b), where 'Hkin ■= (^/94m) = {Cyl{Jl)) by straightforward 
verification and the (•) denotes the completion with respect to the natural equipped inner product 

on Hkin, 

< [a]\[b] >kin-= oJkinia -b). 

To show how this inner product works, given any two cylindrical functions if/ = [iAff],tA' = 
[ilf'a'] 6 Cyl(Jl), the inner product between them is expressed as 

<mm>kin:= f (K,M(P:.Mdfia", (23) 

for any groupoid a" containing both a and a'. The measure d/^a on Xa is defined by the pull 
back of the product Haar measure d/j.^" on the product group S U (2)^" via the identification 
bijection between and 5 ?7(2)^", where Na is number of maximal analytic edges generating 
a. In addition, a nice result shows that given such a family of measures {fia]as£, a probability 
measure /u is uniquely well-defined on the quantum configuration space Jl [fTTl . such that the 
kinematical Hilbert space 'Hkin coincides with the collection of the square-integrable functions 
with respect to the measure // on the quantum configuration space, i.e. Hkin = ^^(^, d/j.), just 
as we expected at the beginning of our construction. 
The representation map nki„ is defined by 

^kin(ci)[b] :=[a-b], V a e and [b] e Hkin- 

Note that nkin{a) is an unbounded operator in general. It is easy to verify that 

nkin{Yf{S)m=m{Y;{S)il,] 
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via Eq. (|22)) . which gives the canonical momentum operator. In the following context, we denote 
the operator nun{Yf{S)) by Pf{S) on "Kfan, and just denote the elements [0-] in l-lkin by ijj for 
simplicity. 

Moreover, since Q.kin '■= 1 is a cyclic vector in "Kyt,,,, the positive linear functional which we 
begin with can be expressed as 

(^kin(o) =< Q.kin\^kin(^)^kin >kin ■ 

Thus the Ashtekar-Isham-Lewandowski state a»yt„, on % is equivalent to a cyclic representation 
CHkin, ^kin, ^kin) for % which is the Ashtekar-Isham-Lewandowski representation for quantum 
holonomy-flux *-algebra of background independent gauge field theory. One thus obtains the 
kinematical representation of loop quantum gravity via the construction of algebraic quantum 
field theory. It is important to note that the Ashtekar-Isham-Lewandowski state is the unique 
state on the quantum holonomy-flux *-algebra % invariant under internal gauge transformations 
and spatial difi'eomorphismfl which are both automorphisms ag and on 31 and can be veri- 
fied that ojki,, ° ag = (i>kw and ojkw ° = cotm- So these gauge transformations are represented 
as unitary transformations on 'Htm, while the cyclic vector Q.kin, representing "no geometry vac- 
uum" state, is the unique state in l-lkm invariant under internal gauge transformations and spatial 
diffeomorphisms. This is a very crucial uniqueness theorem for canonical quantization of gauge 
field theory [|96l : 

Theorem 3.4.1: There exists exactly one Yang-Mills gauge invariant and spatial diffeomor- 
phism invariant state (positive linear functional) on the quantum holonomy-flux * -algebra. In 
other words, there exists a unique Yang-Mills gauge invariant and spatial diffeomorphism in- 
variant cyclic representation for the quantum holonomy-flux *-algebra, which is called Ashtekar- 
Isham-Lewandowski representation. Moreover, this representation is irreducible with respect to 
an exponential version of the quantum holonomy-flux algebra (deflned in M30^ ). which is anal- 
ogous to the Weyl algebra. 

Hence we have finished the construction of kinematical Hilbert space for background inde- 
pendent gauge field theory and represented the quantum holonomy-flux algebra on it. Then 
following the general programme presented in the last subsection, we should impose the con- 
straints as operators on the kinematical Hilbert space since we are dealing with a gauge system. 

3.5 Spin-network Decomposition of Kinematical Hilbert Space 

The kinematical Hilbert space l-lkm for loop quantum gravity has been well defined. In this 
subsection, it will be shown that l-ikm can be decomposed into the orthogonal direct sum of 1- 
dimensional subspaces and find a basis, called spin-network basis, in the Hilbert space, which 
consists of uncountably infinite elements. So the kinematic Hilbert space is non-separable. In 
the following, we will show the decomposition in three steps. 

• Spin-network Decomposition on a Single Edge 

'The proof of this conclusion depends on the compact support property of the smear functions /' (see f9^ for 
detail). 
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Given a graph consisting of only one edge e, which naturally associates with a group 
SU{2) = Xa(e), the elements of Xa(e) are the quantum connections only taking non- 
trivial values on e. Then we consider the decomposition of the Hilbert space Haie) = 
L^(Xa(e),dfia(e)) - L^{S U(2), dfin), which is nothing but the space of square integrable 
functions on the compact group SU{2) with the natural inner product. It is natural 
to define several operators on 'Haie)- First, the so-called configuration operator f(A(e)) 
whose operation on any i/r in a dense domain of L^{S U{2), dfin) is nothing but multipli- 
cation by the function f(A(e)), i.e., 

/(AieMiAie)) := /(A(e))^(A(e)), 

where A(e) 6 SU{2). Second, given any vector^ 6 su(2), it generates left invariant vector 
field L^^' and right invariant vector field R^-^^ on S U (2) by 

L(fV(^(e)) := ^\r=o>ff(A(e)cxpm, 
at 

i?<fV(A(e)) := 4r=o'A(exp(-?^)A(e)), 
at 

for any function i/r e C^(S U(2)). Then one can define the so-called momentum operators 
on the single edge by 

jf' = iL^^-^ and jf^ = iR^^-\ 

where the generators r, e su(2) constitute an orthonormal basis with respect to the 
Killing-Cartan metric. The momentum operators have the well-known commutation re- 
lation of the angular momentum operators in quantum mechanics: 

Third, the Casimir operator on "Kg can be expressed as 

P := = 6'^J^''^J^^\ (24) 

^ J ^ J 

The decomposition of 'He = l}{SU{2),dnH) is provided by the following Peter- Weyl 
Theorem. 

Theorem 3.5.1 (531: 

Given a compact group G, the function space L}{G,diXH) can be decomposed as an or- 
thogonal direct sum of finite dimensional Hilbert spaces, and the matrix elements of the 
equivalence classes of finite dimensional irreducible representations of G form an orthog- 
onal basis in L}{G, dpn)- 

Note that a finite dimensional irreducible representation of G can be regarded as a matrix- 
valued function on G, so the matrix elements are functions on G. Using this theorem, one 
can find the decomposition of the Hilbert space: 

L\SU{2),dnH) = ®jVHj®m, 
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where j, labelling irreducible representations of S U{2), are the half integers, 'Hj denotes 
the carrier space of the j-representation of dimension 2j+l, and "K* is its dual space. The 

basis {ein ® ef } in 'Hj ® 'H* maps a group element g e S U(2) to a. matrix {njnnig)}, where 
m,n = -J, J. Thus the space "Ky ® "H* is spanned by the matrix element functions nim 
of equivalent j-representations. Moreover, the spin-network basis can be defined. 



Proposition 3.5.1 11561 

The system of spin-network functions on ^a(e), consisting of matrix elements {7r/„„} infinite 
dimensional irreducible representations labelled by half-integers [f], satisfies 

p-j^j =/f/+lW' /^V-' =mn^ /'^V-' =nn^ 

where j is called angular momentum quantum number and m,n = magnetic 
quantum number. The normalized functions { -yjl j + lniin}form an orthonormal basis in 
'Ha{e) by the above theorem and 

which is called the spin-network basis on 'Ha{e)- So the Hilbert space on a single edge has 
been decomposed into one dimensional subspaces. 



Note that the system of operators {7 , } forms a complete set of commutable op- 

erators in 'Haie)- There is a cyclic "vacuum state" in the Hilbert space, which is the 
(j = 0)-representation Q^ie) = ^^""^ = 1, representing that there is no geometry on the 
edge. 

• Spin-network Decomposition on Finite Graph 

Given a groupoid a generated by a graph y with A'^ oriented edges and M vertices, 
one can define the configuration operators on the corresponding Hilbert space '74, = 
L\X,, dpj - L2(5 U(2f, Jyu^) by 

/(A(e,)¥(A(ei), A(e;v)) := /(A(e,)¥(A(ei), ...Aie^)). 

The momentum operators 7, ' associated with a edge e connecting a vertex v are defined 
as 

y/""'"^ := (1 (8)... 7;- (8>... (8)1), 
where we set 7,- = if v = s{e) and 7,- = 77^ if v = tie), so J ■ = iXf'"\ Note that 

^ (s v) 

Ji ' only acts nontrivially on the Hilbert space associated with the edge e. Then one 
can define a vertex operator associated with vertex v in analogy with the total angular 
momentum operator via 
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where 

e at V 

Obviously, 'Ha can be firstly decomposed by the representations on each edge e of a as: 

where j := (ji, Ja?) assigns to each edge an irreducible representation of SU(2), in 
the fourth step the Hilbert spaces associated with the edges are allocated to the vertexes 
where these edges meet so that for each vertex v, 

f<;;;'''^^siey-.'K] and 'h;=;^^ ^ ®,^=,^p 

The group of gauge transformations g{v) e SU{2) at each vertex is reducibly represented 
on the Hilbert space 'H^'~f''^ ® 'J^'u!'''^^ in a natural way. So this Hilbert space can be 
decomposed as a direct sum of irreducible representation spaces via Clebsch- Gordon 
decomposition: 

As a result, 'Ha can be further decomposed as: 

'Ha = ®i [^vCez-Hj;,) ,)] = Sj [eK^v-H;;,) /)] = ej [e.'K.j,,] . (25) 

It can also be viewed as the eigenvector space decomposition of the commuting operators 
[Pf (with eigenvalues /(/ + 1)) and [Pf = 6'jJfP. Note that 1 := (Iu-Jm) assigns 
to each vertex(objective) of a an irreducible representation of SU(2). One may also 
enlarge the set of commuting operators to further refine the decomposition of the Hilbert 
space. Note that the subspace of 'Ha with 1 = is Yang-Mills gauge invariant, since the 
representation of gauge transformations is trivial. 

• Spin-network Decomposition of 'Hkin 

Since 'Hkin has the structure 'Hun = ( ^aezHa ), one may consider to construct it as a 
direct sum of 'Ha by canceling some overlapping components. The construction is pre- 
cisely described as a theorem below. 

Theorem 3.5.2: 

Consider assignments j = (ji, jw) to the edges of any groupoid a e £. and assign- 
ments 1 = {Ii,...,Im) to the vertices. The edge representation j is non-trivial on each 
edge, and the vertex representation I is non-trivial at each spuriou^ vertex, unless it is 
the base point of a close analytic loop. Let 'H'^ be the Hilbert space composed by the 

'A vertex v is spurious if it is bivalent and e o e' is itself analytic edge with e, e' meeting at v. 
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subspaces (assigned the above conditions) according to Then IHkm can be 

decomposed as the direct sum of the Hilbert spaces IH'^, i.e., 

"Hkin = ©aexK C. 



Proof: 

Since the representation on each edge is non-trivial, by definition of the inner product, it 
is easy to see that 'H'j, and 'H^, are mutual orthogonal if one of the groupoids a and a' has 
at leat an edge e more than the other due to 

Kindf^e = L 1 • Kndf^e = 

for any j 0. Now consider the case of the spurious vertex. An edge e with j- 
representation in a graph is assigned the Hilbert space Tij <8) 'HJ*. Inserting a vertex v 
into the edge, one obtains two edges ei and e2 split by v both with j-representations, 
which belong to a different graph. By the decomposition of the corresponding Hilbert 
space, 

J J J J J ' J 

the subspace for all / are orthogonal to the space 'HJ <8> 'TYJ*, while the subspace for 
I = coincides with 9ij <8) "Hj* since 'K/^g = C and A(e) = A{ei)A{e2). This completes 
the proof. 

Since there are uncountably many graphs on S, the kinematical Hilbert 'Hkin is non- 
separable. We denote the spin-network basis in 'Hkin by H^, s = (yC*), j^, m^, and 
vacuum Qjt,„ = Eq = 1 , where 

n, := Y\ V2j; + 1<„, o; ^0), 

eeE(y(s)) 

which form a orthonormal basis with the relation < n^in,' >km= ^.v.v'- And Cyly{M) c 
Cyl{Ji) denotes the linear span of the spin network functions for y{s) = y. 

The spin-network basis can be used to construct the so-called spin network representation 
of loop quantum gravity. 

Definition 3.5.1: The spin-network representation is a vector space 'H of complex valued func- 
tions 

^: 5 ^ C; 5 ^ ¥(5), 

where S is the set of the labels s for the spin network states. 'H is equipped with the scalar 
product 

< >:= 
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between square summable functions. 

The relation between the Hilbert spaces 'H and l-lkm is clarified by the following proposition 

Proposition 3.5.2: 

The spin-network transformation 

T : nHun ^ -K; »F ^ ¥(5) :=< n,, ¥ 
is a unitary transformation with inverse 

r-^^ = ^ *i'(^)n,. 

Thus the connection representation and the spin-network representation are "Fourier trans- 
forms" of each other, where the role of the kernel of the integral is played by the spin-network 
basis. Note that, in the gauge invariant Hilbert space of loop quantum gravity which we will 
define later, the Fourier transform with respect to the gauge invariant spin-network basis is the 
so-called loop transform, which leads to the unitary equivalent loop representation of the theory 

imcaimii. 

To conclude this subsection, we show the explicit representation of elementary observables 
on the kinematical Hilbert space ^un- The action of canonical momentum operator Pf(S) on 
differentiable cylindrical functions ifry e Cyly{^) can be expressed as 

Pf(S)if^yi{Aie)UEir)) = I J] f'(v)[J]K{S, e)jy''-^]il^yi{A(e)UEiy)) 

vey(7)n5 e at V 

= I J] /WLTlf^f - 7<f,;Vr({^(^)}.ezr(r))> (26) 

veV{y)nS 

where 

f(S,v) ^ fieuv) f{eu,v) 
•'iiu) - •'i ••• •'i ' 

j(5,v) _ 7<^-..v) ^ _ _ ^ jt.„,,,v)^ ^27) 

for the edges ei, e„ lying above S and eu+d lying below S . And it was proved that 

the operator Pf(S) is essentially self-adjoint on 'Hkin H 14611 . On the other hand, it is obvious to 
construct configuration operators by spin-network functions: 

U,ll/yi{Aie)}eeEiy)) ■= Rs{{Aie)}eeE{r{s)))'ffyi{A{e)}eeE(y))- 

Note that H, may change the graph, i.e., tl/. Cyly(Jl) Cylyuy(s)(Jl). So far, the elementary 
operators of quantum kinematics have been well defined on l-lkin- 
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3.6 Quantum Riemannian Geometry 



The well-established quantum kinematics of loop quantum gravity is now in the same status 
as Riemannian geometry before the appearance of general relativity and Einstein's equation, 
giving general relativity mathematical foundation and offering living place to the Einstein equa- 
tion. Instead of classical geometric quantities, such as length, area, volume etc., the quantities 
in quantum geometry are operators on the kinematical Hilbert space "Hkm, and their spectrum 
serve as the possible values of the quantities in measurements. So far, the kinematical quan- 
tum geometric operators constructed properly in loop quantum gravity include length operator 
[I141L area operator [125] [21], two different volume operators [18] [125] [22], Q operator [IIOOII , 
etc.. Recently, a consistency check was proposed for the different regularizations of the volume 
operator [[TTlirTSl . We thus will only introduce the volume operator defined by Ashtekar and 
Lewandowski ["22^, which is shown to be correct in the consistency check. 

First, we define the area operator with respect to a 2-surface S by the elementary operators. 
Given a closed 2-surface or a surface S with boundary, we can divide it into a large number 
of small area cells S /. Taking account of the classical expression of an area, we set the area of 
the 2-surface to be the limit of the Riemannian sum 

N . 

As := lim[As];v = lim /c^S V Jp^iS j)PjiS j)6'J . 

N~>ca N—>ca i—i \ 

7=1 

Then one can unambiguously obtain a quantum area operator from the canonical momentum 
operators Pi{S) smeared by constant functions. Given a cylindrical function \py 6 Cyly{JK) 
which has second order derivatives, the action of the area operator on i^y is defined in the limit 
by requiring that each area cell contains at most only one intersecting point v of the graph y and 
5 as 

^ \ : 

A5(Ar := lim[A5]/viAy = lim a-^S V JPi(S i)Pj{S ,)6'j iffy. 

N^oo N~>ca Z— i > 

7=1 

The regulator is easy to remove, since the result of the operation of the operator P,(5 /) does 
not change when S j shrinks to a point. Since the refinement of the partition does not affect the 
result of action of \_As]n on ij/y, the limit area operator As, which is shown to be self-adjoint 
|[2n . is well defined on 'Hkin and takes the explicit expression as: 

As^y = Anptl X V(^f-^')(^'-^S>^- 

where f.^^^ and f"^^^ have been defined in Eq. (|27]) . It turns out that the finite linear combinations 
of spin-network basis in l-lkm diagonalizes A5 with eigenvalues given by finite sums, 

as = AnH^l 2 V^^f + D + 2ff\jf + 1) - + 1), (28) 

V 

where f"\ /'^^ and are arbitrary half-integers subject to the standard condition 

jiu.,) g _ |^-(«) _ ^-^1 + 1^ ^ j(u) + (29) 
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Hence the spectrum of the area operator is fundamentally pure discrete, while its continuum 
approximation becomes excellent exponentially rapidly for large eigenvalues. However, in fun- 
damental level, the area is discrete and so is the quantum geometry. One can see that the 
eigenvalue of As does not vanish even in the case where only one edge intersects the surface at 
a single point, whence the quantum geometry is distributional. 

The form of Ashtekar and Lewandowski's volume operator was introduced for the first time 
in IfTSl , and its detailed properties are discussed in [|22| . Given a region R with a fixed coordinate 
system {x^'}a=i,2,3 it, one can introduce a partition of R in the following way. Divide R into 
small volume cells C such that, each cell C is a cube with coordinate volume less than e and 
two different cells only share the points on their boundaries. In each cell C, we introduce three 
2-surfaces s = (5\5^,5^) such that x" is constant on the surface S". We denote this partition 
(C, s) as Pg. Then the volume of the region R can be expressed classically as 

c 

where 

qc, = ^^e'i'n Pt{S")Pj{S')Pu{S% 

3 ! —ahc ■' 

This motivates us to define the volume operator by naively changing PiiS") to PiiS"): 

= liniV.^, 

C 

qc. = ^-^e^\JtiS")Pj{S')P,{S% 

Note that, given any cylindrical function i/fy e Cyly{Ji), we require the vertexes of the graph 
7 to be at the intersecting points of the triples of 2-surfaces s = (5 \ 5 ^,5^) in corresponding 
cells. Thus the limit operator will trivially exist due to the same reason in the case of the area 
operator. However, the volume operator defined here depends on the choice of orientations for 
the triples of surfaces s = (S^,S^,S^), or essentially, the choice of coordinate systems. So it 
is not uniquely defined. Since, for all choice of s = (S\S^,S^), the resulting operators have 
correct semi-classical limit, one settles up the problem by averaging different operators labelled 
by different s [|22|. The process of averaging removes the freedom in defining the volume 
operator up to an overall constant kq. The resulting self-adjoint operator acts on any cylindrical 
function i/zy e Cyly{S\) as 

Vr ifjy = Kq ^ ^|\qv,y\ ^y, 
veV(a) 



where 



e,e',e" at v 
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here e{e,e',e") = sgn{€ahc^''e'^e"'')\v with e" as the tangent vector of edge e and eabc as the 
orientation of Z. The only unsatisfactory point in the present volume operator is the choice 
ambiguity of Kq. However, fortunately, the most recent discussion shows that the overall un- 
determined constant kq can be fixed to be V6 by the consistency check between the volume 
operator and the triad operator [iTTlirTSll . 
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4 Implementation of Quantum Constraints 



After constructing the kinematical Hilbert space 'Hun of loop quantum gravity, one should im- 
plement the constraints on it to obtain the physical Hilbert space which encodes the complete 
information of quantum dynamics of general relativity, since the Hamiltonian of general rela- 
tivity is a linear combination of the constraints. Recalling the constraints ([7]) in the Hamiltonian 
formalism and the Poission algebra ^ among them, the subalgebra generated by the Gauss 
constraints ^(A) forms a Lie algebra and a 2-sided ideal in the constraints algebra. So in this 
section, we first solve the Gaussian constraints independently of the other two kinds of con- 
straints and find the solution space 'H*^, which is constituted by internal gauge invariant quan- 
tum states. Then, although the subalgebra generated by the diffeomorphism constraints is not 
an ideal in the constraint algebra, we still would like to solve them independently of the scalar 
constraints for technical convenience. After that, the quantum operator corresponding to the 
Hamiltonian constraint(scalar constraint) is defined on the kinematical Hilbert space, and we 
will also discuss an alterative for the implementation of the scalar constraint, which is called 
the master constraint programme by modifying the classical constraint algebra. 

4.1 Solutions of Quantum Gaussian Constraint 

Recall that the classical expression of Gauss constraints reads 

^(A) = r dhA'DaP"! = - r dhP'IDaA' = -P(DA), 

where D^A' = daA' + .^A^lsl^. As the situation of triad flux, the Gauss constraints can be 

defined as cylindrically consistent vector fields Y^a on which act on any cylindrical function 

fy e Cylym by 

YdA o fyme)}eeEiy)) := {-P{DA),fy({A(e)leE(y))}. 

Then the Gauss constraint operator can be defined in analogy with the momentum operator, 
which acts on fy as: 

0iA)fy({A(e)}eeE(y)) ■= itiYDAfy({A(e)}eeE(j)) 

= ^ ^ [A'(v)jnf({Me)}eeEM), 

veViy) 

which is the generator of internal gauge transformations on Cyly(Jl). The kernel of the opera- 
tor is easily obtained in terms of the spin-network decomposition, which is the internal gauge 
invariant Hilbert space: 

"f^^ = ®ff,j'^a,j,l=() ® C. 

One then naturally gets the gauge invariant spin-network basis T^, s = (7(5), j s, i*) in "K*^ via a 
group averaging technique at each vertex HI 2611 ll25l Il32l (we will call spin-network state in the 
following context): 

veV(y) eeEiy) 
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assigning a non-trivial spin representation j on each edge and a invariant tensor i (intertwiner) on 
each vertex. We denote the vector space of finite linear combinations of vacuum state and gauge 
invariant spin-network states Cyl{^IQ), which is dense in 'K'^. And Cylyi^lQ) c Cyl(Jl/Q) 
denotes the linear span of the gauge invariant spin network functions for 'y{s) = y. All 
Yang-Mills gauge invariant operators are well defined on 'K'^. However, the condition of acting 
on gauge invariant states often changes the structure of the spectrum of quantum geometric 
operators. For the area operator, the spectrum depends on certain global properties of the surface 
5 (see [1201 [ETI for details). For the volume operators, a non-zero spectrum arises from at least 
4-valent vertices. 

4.2 Solutions of Quantum Diffeomorphism Constraint 

Unlike the strategy in solving Gaussian constraint, one cannot define an operator for the quan- 
tum diffeomorphism constraint as the infinitesimal generator of finite diffeomorphism trans- 
formations (unitary operators since the measure is diffeomorphism invariant) represented on 
"Hkin- The representation of finite diffeomorphisms is a family of unitary operators l]^ acting on 
cylindrical functions \py by 

IJ^ilfy := (Ai^oy, (30) 

for any spatial diffeomorphism on 2. An 1 -parameter subgroup i^, in the group of spatial 
diffeomorphisms is then represented as an 1 -parameter unitary group (7^, on "K/t,,,. However, 
f/y, is not weakly continuous, since the subspaces 'H'^^^^ and 'H'^^^^^^^ are orthogonal to each 
other no matter how small the parameter t is. So one always has 

I < ^Ayl^yJ^r - < ^y\^y >kin I =< ^yl^f^y >km* 0, (31) 

even in the limit when t goes to zero. Therefore, the infinitesimal generator of U^^ does not 
exist. In the strategy to solve the diffeomorphism constraint, due to the Lie algebra structure of 
diffeomorphism constraints subalgebra, the so-called group averaging technique is employed. 
We now outline the procedure. First, given a colored graph (a graph y and a cylindrical function 
living on it), one can define the group of graph symmetries GSy by 

GSy:=Diffy/TDiffy, 

where Diffy is the group of all diffeomorphisms preserving the colored y, and TDiffy is the 
group of diffeomorphisms which trivially acts on y. We define a projection map by averaging 
with respect to GSy to obtain the subspace in Cyly which is invariant under the transformation 
of GSy-. 

for all cylindrical functions ij/y e *T^'a(yy where Uy is the number of the finite elements of GSy. 
Second, we average with respect to all remaining diffeomorphisms which move the graph y. 
For each cylindrical function ifjy 6 Cyly(Jl/Q), there is an element T](if/y) associated to it in the 
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algebraic dual space Cyl* of Cyl{J{IQ), which acts on any cylindrical function (py' e Cylyi^lQ) 
as: 



^eDifmiDiffy 

It is well defined since, for any given graph y', only finite terms are non-zero in the summation. 
It is easy to verify that ri{i]/y) is invariant under the group action of Diff(L), since 



Thus we have defined a rigging map rj : Cyli^lRIQ) Cyl'^^j^ which maps every cylindrical 
function to a diffeomorphism invariant one, where Cyl'^-^^j- is spanned by vacuum state To = 1 
and rigged spin-network functions T[s] = {ri{Ts)}, [s] = ([7],j,i) associated with difFeomor- 
phism classes [y] of graphs y. Moreover a Hermitian inner product can be defined on Cyl'^.^^ 
by the natural action of the algebraic functional: 



The diffeomorphism invariant Hilbert space 9^Diff is defined by the completion of Cyl'^-^j with 
respect to the above inner product < | >Diff - The difFeomorphism invariant spin-network func- 
tions T[s] form an orthonormal basis in 'Hoi/f- Finally, we have obtained the general solutions 
invariant under both Yang-Mills gauge transformations and spatial diffeomorphisms. 

In general relativity, the problem of observables is a subtle issue due to the diffeomorphism 
invariance [I116IIII119IIII120II . Now we discuss the operators as diffeomorphism invariant ob- 
servables on ^Diff- We call an operator O e Jli'Hkm) a strong observable if and only if 
tj'^du^ = 0, ^ if e Diff(L). We call it a weak observable if and only if O leaves "Hz),// 
invariant. Then it is easy to see that a strong observable O must be a weak one. One notices that 
a strong observable O can first be defined on 'Hoiff by its dual operator O* as 






'7(<Ay)[f^¥>07'] = ^('Ay)[0y']- 



< r](ifry)\r](^y) >Diff= r](ij/y)[(py']. 



((9*Oo,y/)[(A] := ^DiffiOif^l 



then one gets 





= < T](d^l(fyM(f>y) >Diff . 
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Note that the Hilbert space "Kc/y is still non-separable if one considers the C" difFeomorphisms 
with n > 0. However, if one extends the difFeomorphisms to be semi-analytic diffeomotphisms, 
i.e. homomorphisms that are analytic difFeomorphisms up to finite isolated points (which can 
be viewed as an extension of the classical concept to the quantum case), the Hilbert space 'Hoiff 
would be separable [|6l[|20l|. 

4.3 Hamiltonian Constraint Operator 

In the following, we consider the issue of scalar constraint in loop quantum gravity. One may 
first construct a Hamiltonian constraint (scalar constraint) operator in l-ikin or 'Hoi/f, then at- 
tempt to find the physical Hilbert space l-lphys by solving the quantum Hamiltonian constraint. 
However, difficulties arise here due to the special role played by the scalar constraints in the 
constraint algebra Q. First, the scalar constraints do not form a Lie subalgebra. Hence the 
strategy of group averaging cannot be used directly on T-lkm for them. Second, modulo the 
Gaussian constraint, there is still a structure function in the Poisson bracket between two scalar 
constraints: 

{^(A^), Sm} = -"ViiNdbM - MdhN)q°\ (32) 

which raises the danger of quantum anomalyies in quantization. Moreover, the diffeomorphism 
constraints do not form an ideal in the quotient constraint algebra modulo the Gaussian con- 
straints. This fact results in that the scalar constraint operator cannot be well defined on 'Hoi/f, 
as it does not commute with the diffeomorphism transformations t/^. Thus the previous con- 
struction of 'Hoiff does not appear very useful for the final construction of 'Hpkys, which is our 
final goal. However, one may still first try to construct a Hamiltonian constraint operator in T-lkm 
for technical convenience. 

We recall the classical expression of Hamiltonian constraint: 

^(A^) := rd\N^^[e\Fi,-2{l+fi')KlKi,] 
^ Jz y I det q\ 

= SEiN)-2{l+/3^)T{N). (33) 

The main idea of the construction is to first express S{N) in terms of the combination of Poisson 
brackets between the variables which have been represented as operators on 'Hkin, then replace 
the Poisson brackets by canonical commutators between the operators. We will use the volume 
functional for a region 7? c E and the extrinsic curvature functional defined by: 

K := Kp^^d'xPlKl 

A key trick here is to consider the following classical identity of the co-triad e'J^x) [I134II : 

(kB)^ P'^Pl 2 

2 ^Jdet^ Kp 
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where Vr is the volume functional for a neighborhood R containing x. And the expression of 
the extrinsic curvature 1-form K'^ix): 

K(x) = ^{A^,(x),K}. 

Note that K can be expressed by a Poisson bracket between the constant- smeared Euclidean 
Hamiltonian constraint and the total volume of the space 2: 

K = I5'\Se{\),V^]. (34) 
Thus one can obtain the equivalent classical expressions of Se(N) and T(N) as: 

= -4z r d'xNixWTr{Fat(x){A,ix),VR^), 



T(N) = J^d'xN-,^KlKi 



pa ph 

' V „ 



VI det^l 

J^d'xNixWTriiAaix), K}{At{x), K}{A,(x), VrJ), 



where = A^t;, Fat = F'^^^i, Tr represents the trace of the Lie algebra matrix, and c Z 
denotes an arbitrary neighborhood of x 6 E. In order to quantize the Hamiltonian constraint 
as a well-defined operator on 'Hkm, one has to express the classical formula of S{N) in terms 
of holonomies A(e) and other variables with clear quantum analogs. As a first attempt 013411 . 
this can be realized by introducing a triangulation T{e), where the parameter e describes how 
fine the triangulation is, and the triangulation will fill out the spatial manifold E when e 
0. Given a tetrahedron A e T{e), we use {5,(A)},=i 2,3 to denote the three outgoing oriented 
segments in A with a common beginning point v(A) = s{Si{A)), and use a,/A) to denote the 
arc connecting the end points of Si{A) and Sj{A). Then several loops q;,/A) are formed by 
aij{A) := Si{A) o a;y(A) o Sj(Ay^. Thus we have the identities: 

{[ AJKAXVr^J = -Aismr'{AiSiiA)),VR^J + oie), 

J.s,(A) 

Aast(A), K} = -Aisi{A)r'{A(si(A)X K} + o(e). 



J.v,(A) 



/ 



Fab(x) = ^A(aij(A)r'-U{aij(A)) + o(e\ 



where Pij is the plane with boundary atj. Note that the above identities are constructed by taking 
account of internal gauge invariance of the final formula of Hamiltonian constraint operator. So 
we have the regularized expression of S{N) by the Riemannian sum 013411 : 

^ Aer(e) 
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TriA(aij(A)r'A(st(A)y'{A(st(A)),VR^.,,^}), 



T\N) 




S\N) 



Tr(A(5,(A))-' {A(5,(A)), i^}A(5/A))-' {A(5/A)), K} x 

AiskiA)r'{A{sk{A)\VR^.J\ 

S%{N)-2{\ +/3^)T'iN), 



(35) 



such that lim^^Q S%N) = S(N). It is clear that the above regulated formula of S(N) is invari- 
ant under internal gauge transformations. Since all constituents in the expression have clear 
quantum analogs, one can quantize the regulated Hamiltonian constraint as an operator on 'Hkm 
(or "K*^) by replacing them by the corresponding operators and Poisson brackets by canonical 
commutators, i.e.. 



Removing the regulator by 6 ^ 0, it turns out that one can obtain a well-defined limit operator 
on l-lkm (or "K*^) with respect to a natural operator topology. 

Now we begin to construct the Hamiltonian constraint operator in analogy with the classical 
expression (|57l) . All we should do is define the corresponding regulated operators on different 
'K^ separately, then remove the regulator e so that the limit operator is defined on 'Hkin (or 'H'^) 
cylindrically consistently. In the following, given a vertex and three edges intersecting at the 
vertex in a graph y oi ij/y 6 Cylj{Jl/ff), we construct one triangulation of the neighborhood 
of the vertex adapted to the three edges. Then we average with respect to the triples of edges 
meeting at the given vertex. Precisely speaking, one can make the triangulations T(e) with the 
following properties HI 3411 [| 14611 . 

• The chosen triple of edges in the graph y is embedded in a T(e) for all e, so that the vertex 
V of y where the three edges meet coincides with a vertex v(A) in T{e). 

• For every triple of segments (ei, 62, £3) of y such that v = s{e\) = s(e2) = ^(^3), there is a 
tetrahedra A 6 T{e) such that v = v(A) = siSi{A)), and Si{A) c e,, V z = 1, 2, 3. We denote 
such a tetrahedra as A° ^ . 

• For each tetrahedra A^^ one can construct seven additional tetrahedron Af, g^.es' P = 
1, 7, by backward analytic extensions of 5,(A) so that t/gi.e,,^, := U^^^Af, g^.es is a neigh- 
borhood of V. 

• The triangulation must be fine enough so that the neighborhoods U(v) := Ug, g^.es Ue, ,62,^3 (v) 
are disjoint for different vertices v and v' of y. Thus for any open neighborhood Uy of the 
graph 7, there exists a triangulation T(€) such that Uy(,v{y)U(v) Q Uy. 

• The distance between a vertex v(A) and the corresponding arcs a;y(A) is described by the 
parameter e. For any two different e and t', the arcs ai/A*^) and aij{A'^') with respect to 
one vertex v(A) are semi-analytically diffeomorphic with each other. 



A(e)^A(e\ ^ V^, {, }^ 
and K ^ = ^[S'e(.1),Vj.I 



[, ] 



ifi 
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• With the triangulations T(€), the integral over E is replaced by the Riemanian sum: 

I - I ■ 

r = w-^ y]^[ +r ]' 

Ju(v) ^.-^^ Jt/„ „„(v) Jj/(v)-C/„ „„,(v) 



where n(v) is the valence of the vertex v = = 5(^2) = ■^'(^3), and ^(v) = ( 1^^^ ) de- 



^ 3 

notes the binomial coefficient which comes from the averaging with respect to the triples 
of edges meeting at given vertex v. One then observes that 



in the limit 6^0. 
• The triangulations for the regions 

U{V) - Uei,e2,e3iv), 

2 - U,, (36) 

are arbitrary. These regions do not contribute to the construction of the operator, since the 
commutator term [A(5,(A)), Vr^,,^,]!/^^ vanishes for all tetrahedron A in the regions (l36l) . 

Thus we find the regulated expression of Hamiltonian constraint operator with respect to the 
triangulations r(e) as [, 134.1 

veV(y) ^ v(A)=v 

Tr(A(5,(A))-i [A(^,(A)), t]A{sjm-' {A{sj{A)), t] x 
A{su{A)r\Ais,{A)\Vu^M 

veViy) 



for any cylindrical function if/y e CylyiJllQ) is a finite linear combination of spin-network states 
r^. with y{s) = y. 
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By construction, the operation of S'^iN) on any i/zy e Cyly{JilQ) is reduced to a finite 
combination of that of SI with respect to different vertices of y. Hence, for each e > 0, S'^{N) 
is a well-defined internal gauge invariant and diffeomorphism covariant operator on Cyl(Jl/^). 

The last step is to remove the regulator by taking the limit 6^0. However, the action of 
the Hamiltonian constraint operator on i/r^ adds arcs a,j(A) with a | -representation with respect 
to each v(A) of 7^, i.e. the action of S'^(N) on cylindrical functions is graph-changing. Hence 
the operator does not converge with respect to the weak operator topology in 'Hkin when e ^ 0, 
since different 'H^(^) with different graphs y are mutually orthogonal. Thus one has to define a 
weaker operator topology to make the operator limit meaningful. By physical motivation and 
the naturally available Hilbert space 'Hoi/f, the convergence of S'^(N) holds with respect to the 
so-called Uniform Rovelli-Smolin Topology [I124II . where one defines S'^(N) to converge if and 
only if ^Diff['S%N)^] converge for all ^Fd,// e ^y^oiff ^ ^ Cyl(Jl/Q). Since the value 
of ^Diff['S'^(.N)(()] is actually independent of e by the fifth property of the triangulations, the 
sequence converges to a nontrivial result ^Diff[<S'^°{N)(p] with arbitrary fixed 60 > 0. Thus we 
have defined a diffeomorphism covariant, densely defined, closed but non- symmetric operator, 
S{N) = limf;^o S'^{N) = S^^iN), on 'Hyt,,, (or "K^) representing the Hamiltonian constraint. 
Moreover, a dual Hamiltonian constraint operator S"^(N) is also defined on Cyl* depending on 
a specific value of e 

(s'%mm := Y[<s^(AO0], 

for all ^ e Cyl* and (p e Cyl(Jl/Q). For ¥0,// e Cyl*^.^^ c Cyl*, one gets 

iS'(N)^'Diff)[cp] = '¥Diff[S%N)cf>]. 

which is independent of the value of e. 

Several remarks on the Hamiltonian constraint operator are listed in the following. 

• Finiteness of S{N) on 'Htin 

In ordinary quantum field theory, the continuous quantum field is only recovered when 
one lets lattice spacing to approach zero, i.e., takes the continuous cut-off parameter to its 
continuous limit. However, this will produce the well-known infinities in quantum field 
theory and make the Hamiltonian operator ill-defined on the Fock space. So it seems 
surprising that our operator S(N) is still well defined, when one takes the limit e ^ 
with respect to the Uniform Rovelli-Smolin Topology so that the triangulation goes to the 
continuum. The reason behind it is that the cut-off parameter is essentially noneffective 
due to the diffeomorphism invariance of our quantum field theory. This is why there 
is no UV divergence in the background independent quantum gauge field theory with 
diffeomorphism invariance. On the other hand, from a convenient viewpoint, one may 
think the Hamiltonian constraint operator as an operator dually defined on a dense domain 
in 'Hoiff- However, we will see that the dual Hamiltonian constraint operator cannot leave 
^Diff invariant. 

''The Hamiltonian constraint operator depends indeed on the choice of the representation j on the arcs a,j(A), 
which is known as one of the regularization ambiguities in the construction of quantum dynamics. For the simpUc- 
ity of the theory, one often choose the lowest label of representation j - \- 
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• Implementation of Dual Quantum Constraint Algebra 

One important task is to check whether the commutator algebra (quantum constraint al- 
gebra) among the corresponding quantum operators of constraints both physically and 
mathematically coincides with the classical constraint algebra by substituting quantum 
constraint operators to classical constraint functional and commutators to Poisson brack- 
ets. Here the quantum anomaly has to be avoided in the construction of constraint opera- 
tors (see the discussion for Eq. (fT3l) ). First, the subalgebra of the quantum diffeomorphism 
constraint algebra is free of anomaly by construction: 

which coincides with the exponentiated version of the Poisson bracket between two dif- 
feomorphism constraints generating the transformations (f, (p' e Diff(Y). 

Second, the quantum constraint algebra between the dual Hamiltonian constraint opera- 
tor S'{N) and the finite diffeomorphism transformation on diffeomorphism-invariant 
states coincides with the classical Poisson algebra between 'V{}i) and S{M). Given a 
cylindrical function (py associated with a graph y and the triangulations T{6) adapted to 
the graph a, the triangulations T{ip o e) = <p o T{e) are compatible with the graph (p o y. 
Then we have by definition: 

( - ([.S(A^), Up]y^Diff)[<Pr] 

= i[S'(N), U;]^'Diff)[(f>y] 

= J] {N(v)^Diff[S:(Py] - Niip O v)Woiff[S^;:,CP,oy]] 
veViy) 

= [N{v) - N(ip O V)]^Diff['Sl(f>y] 

veViy) 

= {S'{N-ip*N)^Diff)[cpyl (37) 

Thus there is no anomaly. However, Eq. (|37l) also explains why the Hamiltonian constraint 
operator iS{N) cannot leave 'Hoiff invariant. 

Third, we compute the commutator between two Hamiltonian constraint operators. No- 
tice that 

{S{N),S{M)]iPy 
= [M(v)S(N) - N(v)S(M)]SlcPy 

veV(y) 

= Z Z \-M{v)N{v')-N(v)M{v')]S';SlcPy, 

veViy) v'eViy') 

where y' is the graph changed from y by the action of S{N) or S{M), which adds the arcs 
aij(A) on y, T(€) is the triangulation adapted to y and r(e') adapted to y'. Since the newly 
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added vertices by SI is planar, they will never contributes the final result. So one has 

[S(N),S(M)^y 

J] [M{v)N{v') - N(v)M{v')]S::Slcf>y 

v,v'eV{y),vi=v' 

= \ iM{vmv')-N(vw(vi][s::si-s::si,](f>y 

vyeV{y),vi=v' 

= \ [M{v)N(v') - A^(v)M(v')][(f/^„,,. - %^,,WA](f>y, 

vyeV{y),v^v' 

(38) 

where we have used the facts that [Sf,, <S^'] = for v v'and there exists a diffeomorphism 
iPvy such that Sl,,Sl = t/^ , >S^,<S^. Obviously, we have in the Uniform Rovelli-Smolin 
Topology 

([S(N),S(M)]y^oiff = 

for all ^Diff 6 Cyl'^-^^-. As we have seen in classical expression Eq. (|32l) . the Poisson 
bracket of any two Hamiltonian constraints is given by a generator of the diffeomrophism 
transformations. Therefore it is mathematically consistent with the classical expression 
that two Hamiltonian constraint operators commute on diffeomorphism invariant states, 
as it is presented above. However, as it has been discussed in [1721 [|95l , the domain of 
dual Hamiltonian constraint operator can be extended to a slightly larger space (habitat) 
in Cyl*, whose elements are not necessary diffeomorphism invariant. And it turns out that 
the commutator between two Hamiltonian constraint operators continues to vanish on the 
habitat, which seems to be problematic. Fortunately, the quantum operator corresponding 
to the right hand side of classical Poisson bracket (|32l) also annihilates every state in the 
habitat [1721 . so the quantum constraint algebra is consistent at this level. But it is not 
clear that whether the quantum constraint algebra, especially the commutator between 
two Hamiltonian constraint is consistent with the classical one (|32l) on some larger space 
in Cyl* containing more diffeomorphism variant stateS On the other hand, more works 
on the semi-classical analysis are also needed to test the classical limit of Eq.(l38]) and 
commutation relation (|32l ). The way to do it is looking for some proper semi-classical 
states for calculating the classical limit of the operators. But due to the graph-changing 
property of the Hamiltonian constraint operator, the semi-classical analysis for the Hamil- 
tonian constraint operator and the quantum constraint algebra is still an open issue so far. 

• General Regularization Scheme of the Hamiltonian Constraint 

In [|20l[, a general scheme of regulation is introduced for the quantization of the Hamilto- 
nian constraint, and includes Thiemann's regularization we introduced above as a specific 
choice. Such a general regularization can be summarized as follows: first, we assign a 

'"However, some scholars disagree with such an argument involving the habitat and consider the habitat to be 
unphysical and completely irrelevant (see, e.g. Ref. 111511 '). 
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partition of Z into cells □ of arbitrary shape. In every cell of the partition we define edges 
Sj, J = 1, and loops jS,, / = 1, n^, where n^, may be different for different cells. 
We use 6 to represent the scale of the cell □ Then fix an arbitrary chosen representation 
p of SU{2). This structure is called a permissible classical regulator if the regulated 
Hamiltonian constraint expression with respect to this partition has correct limit when 

6^0. 

Second, we assign the diffeomorphism covariant property and let the partition adapted to 
the choice of the graph. That is, given a cylindrical function i^y e Cyly{J{IQ), we make 
the partition sufficiently refined that every vertex v 6 V{y) is contained in exact one cell 
of the partition. And if (7, v) is diffeomorphic to (y', v') then, for every e and e' , the quin- 
tuple (7, V, □, {sj), (j8;)) is diffeomorphic to the quintuple (7', V , □', {s'j), (jSp), where 
□ and □' are the cells in the partitions with respect to 7 and 7' respectively, containing v 
and v' respectively. 

As a result, the Hamiltonian constraint operator in this general regularization scheme is 
expressed as: 

t;(N) = Yj ^^Yj^''''T^^(PlMsl')][p[A(sj)], K]p[A(s-/^^^ 

veViy) '^f^ IIK 

X p{A{s-^')]\p{A{sK)^, 1>f/,d), 

where O'' and are fixed constants independent of the value of 6, the values of them 
are determined such that the above expressions have correct classical limits. After remov- 
ing the regulator e via diffeomorphism invariance the same as we did above, we obtain a 
well-defined diffeomorphism covariant operator on 'Hkin (or 'H^) in the sense of the Uni- 
form Rovelli-Smolin Topology, or dual-define the operator on some suitable domain in 
Cyl* . Note that such a general scheme of construction exhibits that there is a great deal 
of freedom in choosing the regulators, so that there are considerable ambiguities in our 
quantization for seeking a proper quantum dynamics for gravity, which is also an open 
issue today. 



4.4 Master Constraint Programme 

Although the Hamiltonian constraint operator introduced above is densely defined on 'Hun and 
diffeomorphism covariant, there are still several unsettled problems which are listed below. 

• It is unclear whether the commutator between two Hamiltonian constraint operators re- 
produces the classical Poisson bracket between two Hamiltonian constraints. Hence it is 
unclear if the quantum Hamiltonian constraint produces the correct quantum dynamics 
with correct classical limit [17211 [|95l . 
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• The dual Hamiltonian constraint operator does not leave the Hilbert space 'Hmff invariant. 
Thus the inner product structure of 'Hoiff cannot be employed in the construction of 
physical inner product. 

• Classically the collection of Hamiltonian constraints does not form a Lie algebra. So one 
cannot employ group averaging strategy in solving the Hamiltonian constraint quantum 
mechanically, since the strategy depends crucially on group structure. 

One may see that all above issues come from the properties of the constraint algebra at classical 
level. However, if one could construct an alternative classical constraint algebra, giving the 
same constraint phase space, which is a Lie algebra (no structure functions), where the subalge- 
bra of diffeomorphism constraints forms an ideal, then the programme of solving the constraints 
would be in a much better position. Such a constraint Lie algebra was first introduced by Thie- 
mann in [I149L The central idea is to introduce the master constraint: 



where C{x) is the scalar constraint in Eq.®. One then gets the master constraint algebra: 

mA^), ^(A^')} = ^([A^,A^']), 
mN), M} = 0, 
{M, M} = 0. 



The master constraint programme has been well tested in various examples [[63l[|64lll65l 
||66lll67l . In the following, we extend the diffeomorphism transformations such that the Hilbert 
space ^Hoiff is separable. This separability of 'Hoiff and the positivity and the diffeomorphism 
invariance of M will be working together properly and provide us with powerful functional 
analytic tools in the programme to solve the constraint algebra quantum mechanically. The 
regularized version of the master constraint can be expressed as 



Introducing a partition V of the 3-manifold E into cells C, we have an operator i?^ acting on 
any cylindrical function fy e Cyly{^IQ) in 'H'^ as 



via a family of state-dependent triangulations T{e) on S as we did in the last section, where 
Xc{v) is the characteristic function of the cell C(v) containing a vertex v of the graph y, and the 




(39) 





(40) 



expression of tiy^ reads 




+2(1 +0^) 



4V2 




(41) 
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Note that /?,,' is similar to that involved in the regulated Hamiltonian constraint operator in 
the last section, while the only difference is that now the volume operator is replaced by its 
quare-root in Eq.lHT]). Hence the action of H^^ on fy adds arcs a,/A) with 1/2-representation 
with respect to each v(A) of y. Thus, for each 6 > 0, is a Yang-Mills gauge invariant and 
diffeomorphism covariant operator defined on Cyl{^IQ). The family of such operators can give 
a limit operator He densely defined on 'H^ by the uniform Rovelli-SmoUin topology. Then a 
master constraint operator, M, acting on any ^oiff 6 ^y^mff '^^^ defined as [,88J 

(moiff)[fy] := ^^Jim^ ^^c(^c)V.], (42) 

for any fy is a finite linear combination of spin-network function. Note that H^(H^yfy is 
also a finite linear combination of spin-network functions on an extended graph with the same 
skeleton of y, hence the value of (M^Di//)[/y] is finite for any given ^oiff- Thus M^oiff lies in 
the algebraic dual of the space of cylindrical functions. Furthermore, we can show that M leaves 
the diffeomorphism invariant distributions invariant. For any diffeomorphism transformation (p 
on S, 

= .4™.„ E 5»::;;c',(»::;;c';)Vri 

CeP 

= rl^,^o'^^'ff^T,\^c(H'c)'fyl (43) 

^ CeP 

where in the last step, we used the fact that the diffeomorphism transformation (p leaves the 
partition invariant in the limit P ^ cr and relabel (p(C) to be C. So we have the result 

(U^M^DifM] = (M¥z),//)[/r]. (44) 

So given a diffeomorphism invariant spin-network state T^g-^, the result state M.T[s] must be a 
diffeomorphism invariant element in the algebraic dual of Cyl{^IQ), which means that 

then 

y'— >Z;e,e'— *0 ^— ' £ 

CeP [ii] 

where the cylindrical function Ecef (^c)^^*2 ^ finite linear combination of spin-network 
functions on some graphs y ' with the same skeleton of 7(52) up to finite number of arcs. Hence 
fixing the diffeomorphism equivalence class [5], only for spin-networks $2 lies in finite number 
of diffeomorphism equivalence class the left hand side of the last equation is non-zero. So there 



46 



are also only finite number of classes [si] in the right hand side such that C[.vi] is non-zero. As a 
result, Mr[i] is a finite linear combination of diffeomorphism invariant spin-network states and 
lies in the Hilbert space of diffeomorphism invariant states 'Hoiff for any {s\. And M is densely 
defined on 'Hoi/f- 

Given two diffeomorphism invariant spin-network functions Tis^] and T[s2], one can give the 
matrix elements of M as L88J [.89] 

< r[,.]|M|r[,,] >o iff 
= (Mr[,,])[r,,e[.v,]] 

= liin yiJ- y y 

}\H^^{H'^^)^ T sie[si] >Kin 

= y lin, yiJ_ y y 

\iicTs >Kin< Ts\{H^(;,yTsielsi] >Kin 

= y y ^li"^ 

[.V] veV(r(.ve[i])) 

X Tu^WJs^ceUA^ Yj < ^.^l(^v')^7^s,e[,v,] >Jfm, (45) 
.v,ce[.v,(;] 

where Diffy is the set of diffeomorphisms leaving the colored graph y invariant, GSy denotes 
the graph symmetry quotient group Diffy/TDiffy where T Diffy is the diffeomorphism which 
is trivial on the graph 7, and Uy is the number of elements in GSy. Note that we have used 
the resolution of identity trick in the fourth step. Since only a finite number of terms in the 
sum over spin-networks s, cells C e f , and diffeomorphism transformations ip are non-zero 
respectively, we can interchange the sums and the limit. In the fifth step, we take the limit 
C ^ V and split the sum Y^s ii^to Zie[i]' where [s, c] denotes the diffeomorphism equivalent 
class associated with s. Here we also use the fact that, given y{s) and y{s') which are different up 
to a diffeomorphism transformation, there is always a diffeomorphism ip transforming the graph 
associated with HIT, (v e y{s)) to that of H^,T,> (v' e r(s')) with ^(v) = v', hence T^.^^iH^.T.eU]] 
is constant for different e [i']. 

Since the term < TsKH^'^Ts^eisi] >Km is independent of the parameter €', one can see 
that by fixing a arbitrary family of state-dependent triangulations T{€'), 



y < TsKH^ )^7"sie[ii] >Kin 
.ve[,v] 

= ^ < )^7".vie[ii] >Kin 

= ^ < -^v U<fiTs\Tsie[st] >Kin 
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z 



< f^v^^-i[^)^7'.v|7'.vie[.vi] >Kii 



= (46) 

where (p are the difFeomorphism transformations spanning the difFeomorphism equivalent class 
[s] . Note that the kinematical inner product in above sum is non- vanishing if and only if (p(y{s))) 
coincides with the graph obtained from certain skeleton yisi) by the action of (H^, )^ and v 6 
Viifiyis))), i.e., the scale (p~\e') of the difFeomorphism images of the tetrahedrons added by 
the action coincides with the scale of certain tetrahedrons in y(s) and <p~^iv) is a vertex in y{s). 
Then we can express the matrix elements (|83l) as: 

< r[,,]|M|r[,,] >Diff 

[.V] veV(y{se[s])) 
[.V] veV(y{se[s])) 

From Eq.(l85l) and the fact that the master constraint operator M is densely defined on 'Hoiff, 
it is obvious that M is a positive and symmetric operator in 'n.Diff- Therefore, the quadratic 
form 2m associated with M is closable m4il . The closure of Qu is the quadratic form of 

a unique self-adjoint operator M, called the Friedrichs extension of M. We relabel M to be 
M for simplicity. From the construction of M, the qualitative description of the kernel of the 
Hamiltonian constraint operator in Ref. l|136ll can be transcribed to describe the solutions to the 
equation: M*F£,,yj = 0. In particular, the diff"eomorphism invariant cylindrical functions based 
on at most 2-valent graphs are obviously normalizable solutions. In conclusion, there exists a 
positive and self-adjoint operator M on 'Hoi/f corresponding to the master constraint (l75l) . and 
zero is in the point spectrum of M. 

Note that the quantum constraint algebra can be easily checked to be anomaly free, i.e., 

[M,f/;]=0, [M,M]=0. 

which is consistent with the classical master constraint algebra in this sense. As a result, the dif- 
ficulty of the original Hamiltonian constraint algebra can be avoided by introducing the master 
constraint algebra, due to the Lie algebra structure of the latter. Since zero is in the spectrum of 
M II1401 . the further task is to obtain the physical Hilbert space 'KpM.v which is the kernel of the 
master constraint operator with some suitable physical inner product, and the issue of quantum 
anomaly is represented in terms of the size of T-lphys and the existence of semi-classical states. 
Note that we will see in the next section that the master constraint programme can be straight- 
forwardly generalized to include matter fields [fS91 . We list some open problems in the master 
constraint programme for further research. 

• Kernel of Master Constraint Operator 

Since the master constraint operator M is self-adjoint, it is a practical problem to define 
DID of 'Hoiff. 
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where fi(A) is the spectral measure with respect to the master constraint operator M. It 
is expected that we can identify "KJ^q with the physical Hilbert space. However, such a 
prescription is ambiguous in the case that zero is only in the continuous spectrum, loses 
physical information in the case that zero is an embedded eigenvalue and unambiguous 
only if zero is an isolated eigenvalue in which case however the whole machinery of the 
DID is not needed at all because "KJ^q c 'Hoiff and the physical inner product coincide 
with the kinematical (differomorphism invariant) one [63] . There are some improved pre- 
scriptions also presented in [|63l by decomposing the measure with respect to the spectrum 
types before direct integral decomposition, some ambiguities can be canceled by some 
physical criterion, e.g., a complete subalgebra of bounded Dirac observables should be 
represented irreducibly as self-adjoint operators on the physical Hilbert space, and the 
resulting physical Hilbert space should admits a sufficient number of semiclassical states. 
Nonetheless, due to the complicated structure of the master constraint operator, it is diffi- 
cult anyhow to manage the spectrum analysis and direct integral decomposition. On the 
other hand, for the self-adjointness of the master constraint operator and the Lie-algebra 
structure of the constraint algebra, a formal group averaging strategy was introduced in 
II149II as a more concrete way to get the physical Hilbert space. It is realized by a formal 
rigged map rjphy,: 



Uphys • ^y^Diff ^ ^phys 



f^%hys(f):=£^<e^'^'f\.>, 



where e'^' is a one parameter continuous unitary group on 'Hoiff by the self-adjointness 
of M, and ^phys is a subset of the algebraic dual of Cyl'^.j-j-. It is trivial to see that riphysif) 
is invariant under the (dual) transformation of e'^' . Thus a inner product can be formally 
defined between two algebraic functionals rjpkysif) and riphysif) in ^phys via: 

< ^physiDl^physif) >phys '■= Vphysif)U']^ 



= [ r diu(A)d(A)]<m\f'iO)>^.^, 
Jr 



where we have used the spectrum decomposition with respect to the self-adjoint operator 
M, the operator e'^' is represented by multiplication by a number e''^' on each "KJ, and 
the vector valued function f{A) is the spectrum decomposition representation of state 
/ e 'Hoiff- Although we can see from the above argument that the physical inner product 
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is proportional to the inner product in the fiber Hilbert space 'Kf^Q, unfortunately, the 
factor diu(A)6(A) is divergent when yu has pure point part, e.g. zero is in the discrete 
spectrum of M. That is one reason why we claim that the above argument is formal. 

On the other hand, the group averaging strategy and the formal physical inner product 
we just defined has potential relationships with path-integral formulation and spin foam 
models due to the positivity of the master constraint operator M I1149L and hopefully, 
we may obtain the physical transition amplitude from this physical inner product in the 
future. However, the whole technique of group averaging for solving the master constraint 
is still formal so far, and the rigorous calculations for it has not done yet as far as we know. 

• Dirac Observables 

Classically, one can prove that a function O e C°°(A1) is a weak observable with respect 
to the scalar constraint if and only if 

{O,{O,M}}\m = 0. 

We define (9 to be a strong observable with respect to the scalar constraint if and only if 

{O,M}\m = 0, 



and to be a ultra-strong observable if and only if 

{o,sm\M = o. 

In quantum version, an observable (9 is a weak Dirac observable if and only if O leaves 
f^phys invariant, while O is now called a strong Dirac observable if and only if O commutes 
with the master constraint operator M. Given a bounded self-adjoint operator O defined 
on 'Hoiff, for instance, a spectral projection of some observables leaving 'Hoifj invariant, 
if the uniform limit exists, the bounded self-adjoint operator defined by group averaging 

I r-T 

[O] := lim — dtU{t)-'dU{t) 

T->ca Zl J_j 

commutes with M and hence becomes a strong Dirac observable on the physical Hilbert 
space. 

• Testing the Classical Limit of the Master Constraint Operator 

One needs to construct spatial difi'eomorphism invariant semiclassical states to calculate 
the expectation value and fluctuation of the master constraint operator. If the results 
coincide with the classical values up to % corrections, one can go ahead to finish our 
quantization programme with confidence. 
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5 Quantum Matter Field on a Quantum Background 



In ordinary quantum field theory, the quantum field is defined on a smooth background space- 
time. However, it is expected that the smooth structure of a spacetime may break down at 
Planck scale, so the present treatment of quantum field theory is valid only in a semiclassical 
sense. Thus we would like to modify the formulation of present quantum field theory to make 
it compatible with the quantum theory of gravity(spacetime) which we already established in 
previous sections so as to explore the behavior of the quantum matter field under Planck scale 
and at extremely strong gravitational fields, e.g. inside the black hole or at the early age of the 
universe. 

In the following, an alternative quantization of scalar field will be introduced, the advantage 
of such a quantization scheme is that the quantum scalar field doesn't depend on the background. 
We will also see that the quantization technique for the previous Hamiltonian constraint can be 
generalized to quantize the Hamiltonian of matter fields coupled to gravity. Then it is shown 
that an operator corresponding to the Hamiltonian of the scalar field can be well defined on 
the coupled diff"eomorphism invariant Hilbert space. It is even positive and self-adjoint without 
any divergence. Thus quantum gravity acts exactly as a natural regulator for the quantum scalar 
field in the polymer representation. Moreover, to study the whole dynamical system of the scalar 
field coupled to gravity, a Hamiltonian constraint operator is defined in the coupled kinematical 
Hilbert space. The contribution of the scalar field to the Hamiltonian constraint can be promoted 
to a positive self-adjoint operator. To avoid possible quantum anomalies and find the physical 
Hilbert space, we will also introduce the master constraint programme for the coupled system. 
A self-adjoint master constraint operator is obtained in the diffeomorphism invariant Hilbert 
space, which assures the feasibility of the programme. 



5.1 Polymer-like Representation of a Scalar Field 

We begin with the total Hamiltonian of the gravity coupled with a massless real scalar field 
which is a linear combination of constraints: 

not = A'Gi + N"Ca + NC, 

where A', N" and are Lagrange multipliers, and the three constraints in the Hamiltonian are 
expressed as If30l[l87]| : 

Gi = D,F; := d.P'! + e,/A',Pl (48) 

Ca = P'^Fl,-A'^Gi + ndacf>, (49) 

C = f^ F;P%e-\Fl, - 2(1 +fi')KlKi,] 

+ -j^V^^5'ip"iP){dam<P + ^n\ (50) 
VI det q\ 2 Ium 

here the real number um is the coupling constant, and ^denotes the momentum conjugate to (p: 

d£ aM 



o(p Jy 
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Thus one has the elementary Poisson brackets 



{A[{x),P){y)] = 6l6)6{x,y), 
[(l){x),n{y)] = 5{x,y). 



Note that the second term of the Hamiltonian constraint (l50l) is just the Hamiltonian of the real 
scalar field. 

Then we look for the background independent representation for the real scalar field coupled 
to gravity, following the polymer representation of the scalar field ||27l . The classical config- 
uration space, 1/, consists of all real-valued smooth functions (p on 2. Given a set of a finite 
number of points X = {^i, ...,xn} inZ, a equivalence relation can be defined by: given two scalar 
field 01, 02 6 1^, 4>i ~ 02 if and only if exp[z/l;0i(A:;)] = exp[//ly02(^7)] for all Xi 6 X and all real 
number Aj. Hence we obtain a bijection between 1// ~ and Rx, which is A'^ copies of the Bohr 
compactification of R [132]. Since one can define a projective family with respect to the set of 
point (graph for scalar field), thus a projective limit 'ZY, which is a compact topological space, 
is obtained as the quantum configuration space of scalar field. Next, we denote by CylxCU) the 
vector space generated by finite linear combinations of the following functions of 0: 



where A = (Ai,A2, • ■ ■,Aj^) are arbitrary non-zero real numbers assigned at each point. It is 
obvious that Cylxiti) has the structure of a *-algebra. The vector space CyliHA) of all cylindrical 
functions on H is defined by the linear span of the linear span of Tq = 1 and Tx^a- Completing 

CylCU) with respect to the sup norm, one obtains a unital Abelian C*-algebra CylCU). Thus 
one can use the GNS structure to construct its cyclic representations. A preferred positive linear 

functional ojq on Cyl{tl) is defined by 



Thus one obtains the Hilbert space, 'H^^ which is defined by L^('U,d/j.), of square integrable 
functions on a compact topological measure space U with respect to fi. The inner product can 
be expressed explicitly as: 



where the label c := (X, A) are called scalar- network. 

As one might expect, the quantum configuration space 1/ is just the Gel'fand spectrum of 

CyliU). More concretely, for a single point set Xq = {xq}, Cylxoitl) is the space of all almost 
periodic functions on a real line R. The Gel'fand spectrum of the corresponding C*-algebra 



TxA(P) '■= ]~[exp[//l^-0(xy)], 




which defines a diffeomorphism-invariant faithful Borel measure // on 1/ as 




(51) 




(52) 
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CylxoCU) is the Bohr completion R^q of R [l27l . which is a compact topological space such that 

CylxoCU) is the C*-algebra of all continuous functions on R^r^. Since R is densely embedded in 
RjcQ, Rvo can be regarded as a completion of R. 

It is clear from Eq.dSTI) that an orthonomal basis in "K^^f is given by the scalar vacuum 
To = 1 and so-called scalar-network functions Tc((p), where c = {X, A) and A = (A^, A2, ■ ■ •, /In) 
are non-zero real numbers. So the total kinematical Hilbert space "K/t,,, is the direct product of 
the kinematical Hilbert space "K^^ for gravity and the kinematical Hilbert space for real scalar 
field, i.e., 'Hkin ■= "H^l^^ ® "H^^ . Then the spin-scalar-network state 7,,^ = T,{A) ® Tc{(p) 6 
Cyly{s){^l&) ® Cylx(c){'U) = Cyly(s,c) is a gravity- scalar cylindrical function on graph y{s, c) = 
7(5) U X(c). Note that generally X(c) may not coincide with the vertices of the graph yis). It is 
straightforward to see that all of these functions constitutes an orthonormal basis in 'Htm as 

< TAA) ® TAmM) ® T,{(P) 6,A'c . 

Note that none of Tikim 'J^^in '^kin ^ separable Hilbert space. 

Given a pair (xq, Aq), there is an elementary configuration for the scalar field, the so-called 
point holonomy, 

U(xo,Ao) ■= exp[iAo^(xo)]. 

It corresponds to a configuration operator U(X(),Aq), which acts on any cylindrical function 
iA(0) 6 Cylx(c)m by 

U{x^, Ao)H4>) = Uixo, AoM(f>). (53) 

All these operators are unitary. But since the family of operators U(xo,A) fails to be weakly 
continuous in A, there is no field operator ^(x) on "H^^. The momentum functional smeared on 
a 3-dimensional region 7? c 2 is expressed by 

n(R) := I d^xn(x). 
Jr 

The Poisson bracket between the momentum functional and a point holonomy can be easily 
calculated to be 

{TTiR), U(x,A)} = -iAxRix)U(x,A), 

where XRi^) is the characteristic function for the region i?. So the momentum operator is defined 
by the action on scalar network functions Tc={x,a) as 

n(R)T,{cf>) := in{n(R),T,m = n[J]AjX(Xj)micf>). 

XjSX 

Now we can impose the quantum constraints on l-lkin and consider the quantum dynamics. 
First, the Gauss constraint can be solved independently of 'H^-^, since it only involves the 
gravitational field. It is also expected that the diff"eomorphism constraint can be implemented 
by the group averaging strategy in the similar way as in the case of pure gravity. Given a spatial 
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diffeomorphism transformation (p, a unitary transformation tj^ was induced by (p in the Hilbert 
space l-ikin, which is expressed as 

Then the differomorphism invariant spin-scalar-network functions are defined by group averag- 
ing as 

where Diffy is the set of diffeomorphisms leaving the colored graph y invariant, GSj denotes 
the graph symmetry quotient group Dif fy/TDif fy where T Diffy is the set of the diffeomor- 
phisms which is trivial on the graph y, and Uy is the number of elements in GSy. Following the 
standard strategy in quantization of pure gravity, an inner product can be defined on the vector 
space spanned by the diffeomorphism invariant spin- scalar-network functions (and the vacuum 
states for gravity, scalar and both respectively) such that they form an orthonormal basis as: 

< 7"[i,c]|7"[.v',c'] >Diff-= T[s^c][Ts',c'e[s',c']] = ^[i,c],[i',c'] • (55) 

After the completion procedure, we obtain the expected Hilbert space of diffeomorphism in- 
variant states for the scalar field coupled to gravity, which is denoted by 'Hoi/f- 

5.2 Diffeomorphism Invariant Hamiltonian of a Scalar Field 

In the following discussion, we consider the quantum scalar field on a fluctuating background. 
A similar idea was considered in Ref. l|139ll , where a Hamiltonian operator with respect to a U(l) 
group representation of the scalar field is defined on a kinematical Hilbert space 'Htm' of matter 
coupled to gravity. Then an effective Hamiltonian operator of the scalar field can be constructed 
as a quadratic form via 

< ^ matters ^matteri^^ 4^ matter ''kin' 
■= < il^gravim) ® ijJmatter, H^atter ^l^gravim) ® ll/'^^tter >kin', (56) 

where iffgravim) 6 "K^^f presents a semiclassical state of gravity approximating some classical 
spacetime background m where the quantum scalar field lives. Thus the effective Hamiltonian 
operator H'^£g^(m) of scalar field contains also the information of the fluctuating background 
metric. In the light of this idea, we will construct a Hamiltonian operator Skg for scalar field 
in the polymer-like representation. It turns out that this Hamiltonian operator can be defined 
in the Hilbert space 'Hoiff of diffeomorphism invariant states for scalar field coupled to gravity 
without UV-divergence. So the quantum dynamics of the scalar field is obtained in a diffeo- 
morphism invariant way, which is expected in the programme of loop quantum gravity. Thus, 
here an effective Hamiltonian operator of the scalar field could be extracted in 'Hoiff by defin- 
ing < *I'[,„](A, 0), Skg ^[m]{^, (p) >Diff to be its expectation value on diffeomorphism invariant 
states ^{(f)) of the scalar field, where the diffeomorphism invariant semiclassical state ¥[„,](A) 
represents certain fluctuating geometry with spatial diffeomorphism invariance, and the label 
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[m] denotes the classical geometry approximated by *P[„](A). Moreover, the quadratic proper- 
ties of the scalar field Hamiltonian will provide powerful functional analytic tools in the quan- 
tization procedure, such that the self-adjointness of the Hamiltonian operator can be proved by 

a theorem in functional analysis. 

Then the crucial point is to define an operator corresponding to the Hamiltonian functional 
Skg of the scalar field, which can be decomposed into two parts 



KG — 'SkG,<^ + Skg, Kin, 



where 



Skg,, = r d'x-^S^P",P){dam<l>, 

^ Jz y I det q\ 

1 r 1 

SkG.Kik = 7; I d^X——=Tp'. 

2aM Js yidet^l 
We will employ the following identities: 

where rf^'^ denotes the Levi-Civita tensor tensity and Vu^ is the volume of an arbitrary neigh- 
borhood containing the point x. By using the point-splitting strategy, the regulated version 
of the Hamiltonian is obtained as: 

Skg,, = '^^J^£yj^£xxe(x-y)6'^x 
1 ,„ 1 



-.Pt(x)(dam)^=Pj(y)M(y) 



}{A,W, V^^} 



Tl"'\dammrAK{x), Vlj^}{AAx), y^n) X 



*KG,Km - ^ 



— f d^x7r(x) fd^yn(y)x 

r 3 dotjejiu)) r det(<(w)) 
Js (yc/„0-^' Js (yc/:;,) ' 

^—7 r r d^ynCy) x 

j3M77"^^Tr({A„(M), Vy^}{A,(M), Vy^}{A,(M), Vy^}) X 



I 

J jVTT-'^^^TrCIA^Cw), V^}{A,(w), Vn;:}{A/w), ylV^J)x 

Xe(x - y)Xe(u - x)Xe{w - y). 
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where we denote by = AlTi,Xeix - y) the characteristic function of a box containing x with 
scale 6 such that Ume^oAfe(-* -y)l^^ = S{x - y), and Vt/j is the volume of the box. In order 
to quantize the Hamiltonian Skg as a well-defined operator in the polymer-like representation, 
we have to express the classical formula of the Hamiltonian in terms of elementary variables 
with clear quantum analogs by introducing a triangulation T{e) of S, where the parameter e 
describes how fine the triangulation is. The quantity regulated on the triangulation is required 
to have correct limit when 6^0. Given a tetrahedron A G T{e), we use {>$,(A)},=i,2,3 to denote 
the three outgoing oriented segments in A with a common beginning point v(A) = s(s,(A)) and 
use (3;/ A) to denote the arcs connecting the end points of ^((A) and 5/ A). Then several loops 
a,/ A) are formed by aij{A) := 5,(A) o a,/ A) o 5^(A)~^ Thus we have the identities: 



J.s(A) 

for nonzero A, where ^(^(A)) and ^(■^'(A)) denote respectively the beginning and end points of 
segment 5(A) with scale e associated with a tetrahedron A. Regulated on the triangulation, the 
scalar field Hamiltonian reads 



{ dtKs\t)yj^ ] = -A(s(A)r\A(s(A)),V'Jf } + o(e). 




and 





^"'"\u(v(A),Ar'[U(t(si(A)),A) - U(v(A),A)] x 
A 

Tr(TiA(s„(A)r'{A(sM),V'Jf }A(s„(A))-\A(s„(A)),V'Jf }) X 

v(A) v(A) 

^P^^U(v(A'),Ar'[U(t(st(A')),A) - U(v(A'),A)] x 
A 

Tr(TjA(Sp(A')r'{A(Sp(A')),V'J,' }A(s,(A')r'{A(s,(A')),Vli' }), 

v(AO v(A') 




A"er(e) 




A"'er(e) 




Xe(v(A) - v(A'))Xe(viA") - v(A));r.(v(A"') - v(A)). 



(57) 
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Note that the above regularization is explicitly dependent on the parameter A, which will lead to 
a kind of quantization ambiguity of the real scalar field dynamics in polymer- like representation. 
Introducing a partition P of the 3-manifold 2 into cells C, we can smear the essential "square 
roots" of iS^g ^ and ^S^^ ^^^^ in one cell C respectively and promote them as regulated operators 
in T-lkin with respect to triangulations T(e) depending on spin- scalar-network state T,,,c as 

W^'Cj. _ y Xc(v) y J, 

v€V(y{s,c)) ^ ' i'(A)=v 
veV(y(i,c)) ^ i<A)=v 

where ;^fc(v) is the characteristic function of the cell C, and 

^^Iv := ^^(v)e''""Tr(A(5,(A))-i[A(^KA)), 7^]x 
A(^„,(A))-i[A(aA)), 7^]x 

A(^„(A))-i[A(^„(A)), ^]). (59) 

Both operators in (l58l) and their adjoint operators are densely defined on "K/t,,,. We now give 
several remarks on their properties. 

• Removal of regulator e 

It is not difficult to see that the action of the operator W^ j' on a spin- scalar-network 
function T^^c is graph-changing. It adds finite number of vertices with representation A{v) 
at ?(5,(A)) with distance e from the vertex v. Recall that the action of the gravitational 
Hamiltonian constraint operator on a spin network function is also graph-changing. As 
a result, the family of operators l^^ j' also fails to be weakly converged when 6^0. 
However, due to the difi'eomorphism covariant properties of the triangulation, the limit 
operator can be well-defined via the uniform Rovelli-Smolin topology, or equivalently, the 
operator can be dually defined on diffeomorphism invariant states. But the dual operator 
cannot leave 'Hoiff invariant. 

• Quantization ambiguity 

As a main difi'erence of the dynamics in polymer-like representation from that in U(l) 
group representation P138I , a continuous label A appears explicitly in the expression of 
(1581 ). Hence there is an one-parameter quantization ambiguity due to the real scalar field. 
Recall that the construction of gravitational Hamiltonian constraint operator also has a 
similar ambiguity due to the choice of the representations j of the edges added by its 
action. A related quantization ambiguity also appears in the dynamics of loop quantum 
cosmology [[501 . 
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Since our quantum field theory is expected to be diffeomorphism invariant, we would like to 
define the Hamiltonian operator of polymer scalar field in the diffeomorphism invariant Hilbert 
space 'Hoiff- For this purpose we fix the parameter /I to be a non-zero constant at every point. 
Then what we will do is to employ the new quantization strategy developed in Refs. HI 4911 
and [I140II . We first construct a quadratic form in the light of a new inner product defined in 
Ref. l|140|| on the algebraic dual £)* of the space of cylindrical functions which is spanned by 
spin-scalar-networks T^ ^ (where the family of labels s, c includes the vacuum states for gravity, 
scalar and both). Then we prove that the quadratic form is closed. Note that, although the 
calculation employing this inner product is formal, it can lead to a well-defined expression 
of the desired quadratic form Eq. (|65l) . Since an arbitrary element of is of the form ¥ = 
Yus,c '^s,c < TsJ ■ >kim one can formally define an inner product < • |- >^ on D* via 

< ^P, r >^ := < C,,, < TJ • >un I Yj "^'s'^C < Ts'A ■ >kin>* 



s,c s' ,c 

— 1 

^s,c^y c' ^ s,c\T s' ,c' ^kin 



s',c' 

s,c:s',c' 



VNX[J,c])N([j',c']) 



where the Cantor aleph N denotes the cardinal of the set [s,c]. Note that we exchange the 
coefficients on which the complex conjugate was taken in Ref.fHOl, so that the inner product 

< "i'Diff Woiff >* reduces to < >Diff for ^^y ^Po,// ^ '^mf- Completing the 

quotient with respect to the null vectors by this inner product, one gets a Hilbert space "K^f. Our 
purpose is to construct a quadratic form associated to some positive and symmetric operator in 
analogy with the classical expression of ( 1571) . So the quadratic form should first be given in 
a positive and symmetric version. It is then natural to define two quadratic forms on a dense 
subset of ^Diff c "K* as: 

(61) 



for any two states ^oiff and which are finite linear combinations of ^[^ f], where the dual 

limit operator W''^ of either family of W^'f or W'^^^ in (|58l) is naturally defined on diffeomor- 
phism invariant states as 



W^'^^z),//[r.J = \im^Diff[W''''T,J. (62) 
To show that the quadratic forms are well defined, we write 

s,c 
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Then, by using the inner product (|60l) the quadratic forms in (1611) become 

CeP ' s,c 

Ce?' ^ [.v,c] '^'-'-■^'''-'^ i.cskc] 



.^^^ Ua-M' l^,m^,c]) ^^^^^^^^^^ 



(63) 

Note that, since ^oiff is a finite linear combination of the diff"eomorphism invariant spin-scalar- 
network basis, taking account of the operational property of W'^ there are only a finite number 
of terms in the summation Yj[s,c] contributing to (|63] ). Hence we can interchange Yj[s,c] ^^id 
limyj^x Yjcep ill above calculation. Moreover, for a sufficiently fine partition such that each cell 
contains at most one vertex, the sum over cells therefore reduces to finite terms with respect to 
the vertices of y(s, c). So we can interchange 2i,c6[i,c] and limp^^ Y^cep to obtain: 

QKG,<l,(^Diff, ^Diff) 

^'^ ^ '^'^Li, t j; ^^^^^^^^^ vey(r(.,c)) 



QKG,Kin(^Diff, ^Diff) 



Diff> 

16 1 



m\ /-/ vL J/ s,ce[s,c]veV(y(s,c)) 



(64) 



where the limit ^ ^ S has been taken so that C ^ v. Since given 7(5, c) and 7(5', c') which 
are diff'erent up to a difFeomorphism transformation, there is always a difi'eomorphism (p trans- 
forming the graph associated with W'^'^'T^^c (v e 7(5', c)) to that of W'^'^''Ts'^c' (y' e 7(5', c')) with 
^(v) = v', (W^"'^z)i7/)[7'.s,ce[i,c]] is constant for different (s, c) e [s, c], i.e., all the N([5, c]) terms 
in the sum over (s, c) e [s, c] are identical. Hence the final expressions of the two quadratic 
forms can be written as: 
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^ [s;c] veViyis,c)) 
QKG,Kin(^Diff, *5'z),-//) 

(65) 

Note that both quadratic forms in (|65l) have finite resuhs and hence their form domains are 
dense in 'Hoi/f- Moreover, both of them are obviously positive, and the following theorem will 
demonstrate their closedness. 



Theorem 5.2.1: Both Qkg,<p <^nd QKCKin are densely defined, positive and closed quadratic 
forms on 'Hoiff, which are associated uniquely with two positive self-adjoint operators respec- 
tively on 'Hoiff such that 

QKG.<l>^Diff, = < 'VDiff\SKG,<t!\'^'Diff >Diff 

QKG,Km(^ Diff^^'uiff-) - ^ Diffl'S KG.Kinl^Diff Diff • 

Therefore the Hamiltonian operator 

'SkG '■= SkG,(P + S KG, Kin (66) 

is positive and also have a unique self-adjoint extension. 



Proof: We follow the strategy developed in Refs. 1114011 and IfSSlI to prove that both Qkg4 and 
QKG,Kin are closeable and uniquely induce two positive self-adjoint operators Skg.<i, and SKG,Km- 
One can formally define <Skg.(/> and Skckw acting on difl'eomorphism invariant spin-scalar net- 
work functions via: 

'Skg,(/> T[si,ci] 

Then we need to show that both of the above operators are densely defined on the Hilbert 
space 'Hoiff- Given a diffeomorphism invariant spin-scalar network function T^^i^ci], there are 
only a finite number of terms Ti,^ c,][W'^'^'T s.ce[s,c]] which are nonzero in the sum over equivalent 
classes [s, c] in (|65l) . On the other hand, given one spin-scalar-network function Ts,ce[s,c], there 

are also only a finite number of possible ^[.v^.cz] such that the terms T[s2,c2][W'^'^'T^^^ce[s,c]] are 
nonzero. As a result, only a finite number of terms survive in both sums over [s2, C2] in Eqs. 
(1671) and (|68l) . Hence both Skg,(P and Skg.kw are well defined on spin-scalar-network basis. 
Then it follows from Eqs. (l65l) . (|67l) and (l68l) that they are positive and symmetric operators 
densely defined in "Kd,//, whose quadratic forms coincide with Qkg4 and Qkg,kw on their form 
domains. Hence both Qkg,(I> and Qkg.kih have positive closures and uniquely induce self-adjoint 



2_j QkG,^{Tis2,c2\, Tis^^ci\)Tls2,C2\, 
.«2,C2] 

^ QKG,Kin(T[s2.C2]' '^[si,ci])T[s2,C2]- 



(67) 
(68) 



[.«2,C2] 
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(Friedrichs) extensions of Skg,(P and Skckm respectively [|1141. which we denote by Skg,<i> and 
SKG,Kin as well. As a result, the Hamiltonian operator Skg defined by Eq.(l66l) is also positive 
and symmetric. Hence it has a unique self-adjoint (Friedrichs) extension. 
□ 

We notice that, from a different perspective, one can construct the same Hamiltonian oper- 
ator 'Hkg without introducing an inner product on D*. The construction is sketched as follows. 
Using the two well-defined operators W^^'f and as in (|58l) . as well as their adjoint operators 
(H^'f )^ and {W^^^y , one may define two operators on "Kd;// corresponding to the two terms in 
dSj'by 

(SKG,Kin'i'Diff)[T,J = lim vj/ [y 84x^6 ^ C(^.'x:)t7^ ]^ 

(69) 

for any spin- scalar-network T..^. In analogy with the discussion about the master constraint 
operator and Ref.[88J, it can be shown that both above operators leave 'Hoiff invariant and are 
densely defined on 'Hoiff- Moreover, the quadratic forms associated with them coincide with 
the quadratic forms in (l65l) . Thus the Hamiltonian operator Skg '■= ^kg^ + Skg, Kin coincides 
with the one constructed in the quadratic form approach. 

In summary, we have constructed a positive self-adjoint Hamiltonian operator on 'Hoiff for 
the polymer-like scalar field, depending on a chosen parameter A. Thus there is an 1 -parameter 
ambiguity in the construction. However, there is no UV divergence in this quantum Hamiltonian 
without renormalization, since quantum gravity plays the role of a natural regulator for the 
polymer-like scalar field. 

5.3 Hamiltonian Constraint Equation for the Coupled System 

In this section we consider the whole dynamical system of scalar field coupled to gravity. Re- 
call that in perturbative quantum field theory in curved spacetime, the definition of some basic 
physical quantities, such as the expectation value of the energy-momentum, is ambiguous and 
it is challenging difficult to calculate the back-reaction of quantum fields on the background 
spacetime [,157,1 . This is reflected by the fact that the semi-classical Einstein equation, 

Ra/S[g] - ^R[g]gaf^ = K< t,f,[g] >, (70) 

are known to be inconsistent and ambiguous Ii70llll46 l|. One could speculate that the diffi- 
culty is related to the fact that the usual formulation of quantum field theories are background 
dependent. Following this line of thought, if the quantization programme is by construction 
non-perturbative and background independent, it may be possible to solve the problems fun- 
damentally. In loop quantum gravity, there is no assumption of a priori background metric at 
all. The quantum geometry and quantum matter fields are coupled and fluctuating naturally 
with respect to each other on a common manifold. On the other hand, there exists the "time 
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problem" in quantum theory of pure gravity, since all the physical states have to satisfy certain 
version of quantum Wheeler-DeWitt constraint equation. However, the situation could improve 
when matter field is coupled to gravity [|54l|[|122L In the following construction, we impose the 
quantum Hamiltonian constraint on l-ikin, and thus define a quantum Wheeler-DeWitt constraint 
equation for the scalar field coupled to gravity. Then one can gain an insight into the problem 
of time from the coupled equation, and the back-reaction of the quantum scalar field is included 
in the framework of loop quantum gravity. 

We now define an operator in ^iun corresponding to the scalar field part Skg{N) of the total 
Hamiltonian constraint functional, which can be read out from Eqs. (l35l) and (|50l) as 



where 



Skg{N) = SkgAN)+SkgMN), 



SkgAN) = f d\N—^SiptP){dam<P, 

2 Jy. A/ldetd 



Vldet^l 

^ - ( d-'xN —^ 
2Q'm Jy ^1 det q\ 



SkgMN) = ^ { d'xN ^ ' ^ 



In analogy with the regularization and quantization in the previous section, the regulated version 
of quantum Hamiltonian constraint S^^^iN) of scalar field is expressed by taking the limit C 
v: 

veV(y(s,c)) ^ 

where for any v e V{y{s, c)), the operators 

^ ' v(A)=v 
^ ' v(A)=v 

and their adjoints are all densely defined in 'Hun- Hence the family of Hamiltonian constraint 
operators dTTI ) is also densely defined, and the regulator e can be removed via the Uniform 
Rovelli-SmoUin topology, or equivalently the limit operator dually acts on diffeomorphism in- 
variant states as 

= lim^D,-//[<5^G(A^)/]' (72) 

e— ►() 



for any / 6 CyK^IQ) ® CylinA). Similar to the dual of Sgr{N), the operator S'j^(^{N) fails to 
commute with the dual of finite diff"eomorphism transformation operators, unless the smearing 
function N{x) is a constant function over S. In fact, the dual Hamiltonian constraint opera- 
tor smeared by = 1 is just the diffeomorphism invariant Hamiltonian we just defined in the 
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last subsection. From Eq.dTTI). it is not difficult to prove that for positive N{x) the Hamilto- 
nian constraint operator Skg(N) of a scalar field is positive and symmetric in "7-4;,, and hence 
has a unique self-adjoint extension |[S91l . Our construction of Skg(N) is similar to that of the 
Higgs field Hamiltonian constraint in Ref. lll38L However, like the case of Skg, there is a 
one-parameter ambiguity in our construction of <Skg(N) due to the real scalar field, which is 
manifested as the continuous parameter A in the expression of n^^,. in (|59l) . Note that now A is 
not required to be a constant, i.e., its value can be changed from one point to another. Thus the 
total Hamiltonian constraint operator of scalar field coupled to gravity has been obtained as 

S{N) = Sgr(N) + SKGiN). (73) 

Again, there is no UV divergence in this quantum Hamiltonian constraint. Recall that, in stan- 
dard quantum field theory the UV divergence can only be cured by a renormalization procedure, 
in which one has to multiply the Hamiltonian by a suitable power of the regulating parameter 
6. However, now e has naturally disappeared from the expression of (l73l) . So renormalization is 
not needed for the polymer-like scalar field coupled to gravity, since quantum gravity has played 
the role of a natural regulator. This result heightens our confidence that the issue of divergences 
in quantum field theory can be cured in the framework of loop quantum gravity. 

Now we have obtained the desired matter-coupled quantum Hamiltonian constraint equation 

- (S'^G(N)'VDiff)[f] = (SUm'VDiffM. (74) 
Comparing it with the well-known Schordinger equation for a particle, 

ifi—i]/{x, t) = H{x, -ifi—)i]/{x, t), 
ot ox 

where if/{x, t) 6 L^(R, dx) and ? is a parameter labeling time evolution, one may take the view- 
point that the matter field constraint operator S'^^^iN) plays the role of ifij^. Then (p appears as 
the parameter labeling the evolution of the gravitational field state. In the reverse viewpoint, the 
gravitational field would become the parameter labeling the evolution of the quantum matter 
field. Note that such an idea has been successfully applied in a loop quantum cosmology model 
to help us to understand the quantum nature of big bang in the deep Planck regime [[28ll[|29l . 

5.4 Master Constraint for the Coupled System 

Recall that in order to avoid possible quantum anomalies and find the physical Hilbert space 
of quantum gravity, the master constraint programme was first introduced in the last section. 
The central idea is to construct an alternative classical constraint algebra, giving the same con- 
straint phase space, which is a Lie algebra (no structure functions) and where the subalgebra of 
diffeomorphism constraints forms an ideal. A self-adjoint master constraint operator for loop 
quantum gravity is then proposed on 'Hoiff- The master constraint programme can be gener- 
alized to matter fields coupled to gravity in a straightforward way. We now take the massless 
real scalar field to demonstrate the construction of a master constraint operator according to 
the same strategy as we did in the last section. By this approach one not only avoids possible 
quantum anomalies which might appear in the conventional canonical quantization method, but 
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also might give a qualitative description of the physical Hilbert space for the coupled system. 
We introduce the master constraint for the scalar field coupled to gravity as 



M := ^ I d'x—EM=, (75) 



41' 



Vldet^(j» 



where C{x) is the Hamiltonian constraint in (1501) . After solving the Gaussian constraint, one 
gets the master constraint algebra as a Lie algebra: 

mN),^(N')} = ^([A^,^']), 
{^(A^), M} = 0, 

{M, M} = 0, (76) 

where the subalgebra of diffeomorphism constraints forms an ideal. So it is possible to define a 
corresponding master constraint operator on 'Hoiff- In the following, the positivity and the dif- 
feomorphism invariance of M will be working together properly and provide us with powerful 
functional analytic tools in the quantization procedure. 

The regulated version of the master constraint can be expressed via a point- splitting strategy 

as: 

:= o d y d xXe{x-y)—=—=. {11) 



Introducing a partition "P of the 3-manifold E into cells C, we have an operator acting on 
any spin-scalar-network state Ts,c via a family of state-dependent triangulation T(6), 

veV(y{s,c)) ^ ' v(A)=v 



-I- 

veV(y(i,cO) 



16 



where 



= ^^6'^*^Tr(A(a,y(A))-iA(^,(A))-HA(^,(A)), ^]) 

+ (1 + 0") 3 ■^3^^3 e'Hx{,A{smr' [A(5,(A)), 
X A(5XA))-i[A(5y(A)),^^]A(5,(A))-HA(5,(A)), ^]), 
= ^ Z f"""^t>(v,^rHW^/(A)),i)-f/(v,^)] 

v(A)=v 

X Tr(T,-A(5„(A))-i [A(5„(A)), t/J^']A(5„(A))-HA(^„(A)), Vg']), 



^ 6,V 
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X TT(A(si(A)r'[A(s,(A)), V'JpA(s^(A)r'[A(s^(A)), V'Jp 

X A(s„iA)r'[A(s„iA)),rj,'^]). (79) 

Hence the action of on a cylindrical function fy adds analytical arcs a,y(A) with ^ -representation 
and points at t(si(A)) with representation constant A with respect to each vertex v(A) of y. Thus, 
for each e > 0, is a SU (2) gauge invariant and diffeomorphism covariant operator defined 
on Cyl{J{.IQ) ® CyliHA). The limit operator He is densely defined on 'Hnn by the uniform 
Rovelli-Smolin topology. And the same result holds for the adjoint operator (^^^^ 
Then a master constraint operator, M, on "Hd,// can be defined by: 

{M^oiffW,,,^ := lim^ WiH^n.^^l (80) 

C€!P 

Since H^(H^)''Ts ^ is a finite linear combination of spin- scalar-network functions on an graph 
with skeleton y, the value of {M^oiff)[Ts,c] is finite for a given ^oiff that is a finite linear 
combination of T[,^c]- So M^/),// is in the algebraic dual of the space of cylindrical functions. 
Moreover, we can show that it is diff"eomorphism invariant. For any diffeomorphism transfor- 
mation ^, 

(U;W¥oiff)[fy] = Jm^ J oiff[J]\Hl(H^c)'U,fy^ 

^ ^ C€P 

CeP 

= ^J^^^o^«^^^[Z^^c(^c')Vr], (81) 

for any cylindrical function f-y, where in the last step, we used the fact that the diff"eomorphism 
transformation cp leaves the partition invariant in the limit P ^ H and relabel (p(C) to be C. So 
we have the result 

(C/;M^z>,//)[/r] = (MWDiff)[fyl (82) 

So given a diffeomorphism invariant spin-scalar-network state , c], the result state MT^s,c] must 
be a diffeomorphism invariant element in the algebraic dual of Cyl{J{.IQ)^Cyl{'U), which means 
that 



then 

C€!P [ii,ci] 
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where the cylindrical function licey j^c^^D^^S2,c2 is a finite linear combination of spin-scalar- 
network functions on some graphs y ' with the same skeleton of 7(52, C2) up to finite number of 
arcs and vertices. Hence fixing the diff'eomorphism equivalence class [s, c], only for spin-scalar- 
networks S2, C2 lying in a finite number of diff'eomorphism equivalence class on the left hand 
side of the last equation is non-zero. So there are also only finite number of classes [.^i, ci] in 
the right hand side such that C[si,ci] is non-zero. As a result, MT^s^^] is a finite linear combination 
of diff'eomorphism invariant spin-network states so lies in the Hilbert space of diffeomorphism 
invariant states 'Hoiff for any [s, c]. And M is densely defined on 'Hoiff- 

We now compute the matrix elements of M. Given two diff'eomorphism invariant spin- 
scalar-network functions Tisi^ci] and T[S2,C2]'> the matrix element of M is calculated as 

< T[si,ci]\M\Tls2^c2] >Diff 



= (Mr[,,,,,])[r,, 

,Ci€[5i,Ci]J 

■;^T[s2,C2][Hc(HcyTsi,Cl€[suCl]] 



lirn v-— y y 



,C2) 



Cep ^y(s2,C2) ^DifflDiffy^,^,,^) V'^GSy^S2,. 

X < f/^[/^/r^2,C2e[i2,C2]l-^c^-^C-^^'^ii,C2e[ii,ci] ^Kin 

= y yi^ y y 

X < U ^tj ^'T s^^c2e[s2,c2\\H-'c^ s,c >Kin< ^sAi^cy^T^ suci^[suci\ >Kin 

= y y iiim 

[s,c] v€V(y(s,c€[s,c])) 

,c€[s,c]\ <Ts,c\(HyyTs^^cl€[sl,Cl]>Kim (83) 

s,ce[s,c] 

where Diffy is the set of diff'eomorphisms leaving the colored graph y invariant, GSy denotes 
the graph symmetry quotient group Diffy/T Diffy where T Diffy is the difteomorphisms which 
is trivial on the graph y, and Uy is the number of elements in GSy. Note that we have used the 
resolution of identity trick in the fourth step. Since only a finite number of terms in the sum 
over spin-scalar-networks (s, c), cells C e P, and diff'eomorphism transformations (p are non- 
zero respectively, we can interchange the sums and the limit. In the fifth step, we take the 
limit C ^ V and split the sum 2.s,c into 1j[s,c] lii,c€[i,c]' where [s, c] denotes the diffeomorphism 
equivalence class associated with (s, c). Here we also use the fact that, given y(s, c) and y(s', c') 
which are diff'erent up to a diffeomorphism transformation, there is always a diflreomorphism 
(f transforming the graph associated with H^y^^^^^T.^c (v e y(s,c)) to that of (v' e 

y(s', c')) with ^(v) = v', hence Ti,^,c2][H'y^^^^^T,^c^y^c]\ is constant for diff'erent (s, c) e [s, c]. 

Since the term 2.v,re[.s,c] < TsJ(Hf'yTs^^c,e[s,,ci] >Km is independent of the parameter e', one 
can see that by fixing a family of arbitrary state-dependent triangulations T(6'), 

y < TsA(Hy )^7'^i,ci€[ii,ci] >Km 

S,C€[S,C] 
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— ^ \ ^ f^(/)7^i,c|7^ii,cie[ii,ci] ^Kin 



- 7'[„,ci][^J!'gy(^(,,^))r,,c], (84) 

where (f are the diffeomorphism transformations spanning the diffeomorphism equivalence class 
[s,c]. Note that the kinematical inner product in the above sum is non- vanishing if and only 
if (fiyis, c))) coincides with the extended graph obtained from certain skeleton y(si,Ci) by the 
action of {H^,'y and v 6 V((p{y(s, c))), i.e., the scale (p'^ie') of the diffeomorphism images of the 
tetrahedrons added by the action coincides with the scale of certain tetrahedrons in y{s, c) and 
^"'(v) is a vertex in y{s, c). Then we can express the matrix elements (l83l) as: 



[.v,c] veV{y(s,ce[s,c\)) 

2 



ZV~^ 1 ^ ^ ^ 

2_i ■j{H[T[s2,c2l)\Ts,ce[s,c]\{H^T[si )[7'i,ce[.v,c] ] • 



(85) 



[s,c\ veV{y(s,ce[s,c\)) 



From Eq.(l85l) and the result that the master constraint operator M is densely defined on "Hoiff, 
it is obvious that M is a positive and symmetric operator on 'Hoiff- Hence, it is associated with 

a unique self-adjoint operator M, called the Friedrichs extension of M. We relabel M to be 
M for simplicity. In conclusion, there exists a positive and self-adjoint operator M on 'Hoi/f 
corresponding to the master constraint (1751) . It is then possible to obtain the physical Hilbert 
space of the coupled system by the direct integral decomposition of 'Hoi/f with respect to M. 

Note that the quantum constraint algebra can be easily checked to be anomaly free. Eq.(l82l) 
assures that the master constraint operator commutes with finite diffeomorphism transforma- 
tions, i.e., 

[M, t/;] = 0. (86) 

Also it is obvious that the master constraint operator commutes with itself, 

[M, M] = 0. (87) 

So the quantum constraint algebra is precisely consistent with the classical constraint algebra 
([761) in this sense. As a result, the difficulty of the original Hamiltonian constraint algebra can 
be avoided by introducing the master constraint algebra, due to the Lie algebra structure of the 
latter. 
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6 The Semiclassical Limit of Quantum Dynamics 



As shown in previous chapters, both the Hamiltonian constraint operator S(N) and the master 
constraint operator M can be well defined in the framework of loop quantum gravity. However, 
since the Hilbert spaces 'Hkm and 'Noiff, the operators S(N) and M are constructed in such 
ways that are drastically different from usual quantum field theory, one has to check whether 
the constraint operators and the corresponding algebras have correct semiclassical limits with 
respect to suitable semiclassical states. 

6.1 The Construction of Coherent States 

In order to find the proper semiclassical states and check the classical limit of the theory, the 
idea of a non-normalizable coherent state defined by a generalized Laplace operator and its 
heat kernel was introduced for the first time in II26L Recently, kinematical coherent states were 
constructed in two different approaches. One leads to the so-called complexifier coherent states 
proposed by Thiemann et al [ 142 l | HI 43l[|144| 11145 L The other was proposed by Varadarajan 
lfT5l]llfT53ll [fT541 and further developed by Ashtekar etal ifTOlllfBll . 

The complexifier approach is motivated by the coherent state construction for compact Lie 
groups lf86l . One begins with a positive function C (complexifier) on the classical phase space 
and arrives at a "coherent state" i/',,,, which more possibly belongs to the dual space Cyl* rather 
than "Kytm- However, one may consider the so-called "cut-off state" of i/^„, with respect to a 
finite graph as a graph-dependent coherent state in 'Hkin II146II . By construction, the coherent 
state i/ffn is an eigenstate of an annihilation operator coming also from the complexifier C and 
hence has the desired semiclassical properties [I143II 111441 . We now sketch the basic idea of 
its construction. Given the Hilbert space 'H for a dynamical system with constraints and a 
subalgerba of observables S in the space JICH) of linear operators on 'H, the semiclassical 
states with respect to S are defined in Definition 3.1.5. Kinematical coherent states {^PmlmeM 
are semiclassical states which in addition satisfy the annihilation operator property II142IIII146II , 
namely there exists a certain non-self-adjoint operator z = a + iAb with a, b e S and a certain 
squeezing parameter A, such that 

|¥,„=z(m)¥,„. (88) 

Note that Eq.dFSl) implies that the minimal uncertainty relation is saturated for the pair of ele- 
ments (a, b), i.e., 

^M([a - ^™(a)]2) = vF„,([S - Y,„(^)]') = b])\. (89) 

Note also that coherent states are usually required to satisfy the additional peakedness property, 
namely for any m e Al the overlap function | < *F„,, > I is concentrated in a phase volume 
5l*I'm([4. ^])L where § is a configuration operator and p a momentum operator. So the central 
element in the construction is to define a suitable "annihilation operator" z in analogy with the 
simplest case of harmonic oscillator. A powerful tool named as "complexifier" is introduced in 
Ref. [|142 | | to define a meaningful z operator which can give rise to kinematical coherent states 
for a general quantum system. 
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Definition 6.1.1: Given a phase space M = T*Cfor some dynamical system with configuration 
coordinates q and momentum coordinates p, a complexifier, C, is a positive smooth function on 
Ai, such that 

(1) C/h is dimensionless; 

(2) lim||p||^oo '^j^ = CO for some suitable norm on the space of the momentum; 

(3) Certain complex coordinates (z(m),z(m)) ofM can be constructed from C. 

Given a well-defined complexifier C on phase space M, the programme for constructing co- 
herent states associated with C can be carried out as the following. 

• Complex polarization 

The condition (3) in Definition 7.3.1 implies that the complex coordinate z(m) of Ai can 
be constructed via 

oo 

z(m):=y-W,C}(„)(m), (90) 
to 

where the multiple Poisson bracket is inductively defined by {q, C}(0) = q, {q, C}(„) = 
{{q,C}(„-i),C}. One will see that z(m) can be regarded as the classical version of an 
annihilation operator. 

• Defining the annihilation operator 

After the quantization procedure, a Hilbert space "K = L^(C, d/u) with a suitable measure 
dfi on a suitable configuration space C can be constructed. It is reasonable to assume that 
C can be defined as a positive self-adjoint operator C on 'H. Then a corresponding oper- 
ator z can be defined by transforming the Poisson brackets in Eq.(|90l) into commutators, 
i.e., 

^^=Z^74r[^'^](") = ^"'''^^'''' (91) 

which is called as an annihilation operator. 
Constructing coherent states 

Let Sq'(q) be the ^-distribution on C with respect to the measure dfi. Since C is assumed 
to be positive and self-adjoint, the conditions (1) and (2) in Definition 7.3.1 imply that 
g-c/s |g ^ well-defined "smoothening operator". So it is quite possible that the heat kernel 
evolution of the ^-distribution, e~'-'^^6q'{q), is a square integrable function in 'H, which is 
even analytic. Then one may analytically extend the variable q' in e'^l^5q'{q) to complex 
values z{m) and obtain a class of states as 

il/Sq) := Ve-^l''6,.{q)\,^,^,n), (92) 

such that one has 

z if^M ■■= [e~'^''qS,'m,'^z(,n) = Vq'e-^l^6,,{q)\'^z(,n) = z(m) i^'Jq). (93) 
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Hence i/r^ is automatically an eigenstate of the annihilation operator z- So it is natural to 
define coherent states ifrmiq) by normalizing ij/'^iq). 

One may check that all the coherent state properties usually required are likely to be satisfied 
by the above complexifier coherent states i/zmiq) [I146II . As a simple example, in the case of one- 
dimensional harmonic oscillator with Hamiltonian H = + \mo/q^), one may choose the 
complexifier C = p^/{2moj). It is straightforward to check that the coherent state constructed 
by the above procedure coincides with the usual harmonic oscillator coherent state up to a phase 
II146II . So the complexifier coherent state can be considered as a suitable generalization of the 
concept of usual harmonic oscillator coherent state. 

The complexifer approach can be used to construct kinematical coherent states in loop quan- 
tum gravity. Given a suitable complexifier C, for each analytic path e c S one can define 

A^(e):= y-{A(e),C}(„), (94) 

n=0 

where A(e) 6 SU{2) is assigned to e. As the complexifier C is assumed to give rise to a 
positive self-adjoint operator C on the kinematical Hilbert space l-lkin, one further supposes that 
C/hTs = tAsTs, where r is a so-called classicality parameter, {Ts(A)}s form a basis in l-fkm 
and are analytic inAe^. Moreover the (5-distribution on the quantum configuration space Jl 
can be formally expressed as Sa'(A) = 2i Ts{A')Ts(A). Thus by applying Eq.(l92l) one obtains 
coherent states 

(A;c(A) = (e-'^"')6AA)U'^AC = e-''^T,M'')TM). (95) 

However, since there are an uncountably infinite number of terms in the expression (l95l) . the 
norm of lA^cC^) would in general be divergent. So tA^c(A) is generally not an element of 'Hkin 
but rather a distribution on a dense subset of 'Hkin- In order to test the semiclassical limit of 
quantum geometric operators on 'Hkin, one may further consider the "cut-off" state" of iA^c(^) 
with respect to a finite graph y as a graph-dependent coherent state in 'Hkin 1114611 . So the 
key input in the construction is to choose a suitable complexifer. There are vast possibilities 
of choice. For example, a candidate complexifier C is considered in Ref. [148il such that the 
corresponding operator acts on cylindrical functions fy by 

iCmfy = \ij] ^(^)-^1)/r' (96) 

eeEir) 

where is the Casimir operator defined by Eq.(l24l) associated to the edge e, the positive num- 
bers 1(e) satisfying l{e o e') = 1(e) + l{e') and l{e'^) = 1(e) serves as a classicalization parameter. 
Then it can be shown from Eq.(|94l) that A'^{e) is an element of 5L(2, C). So the classical inter- 
pretation of the annihilation operators is simply the generalized complex SU{2) connections. It 
has been shown in Refs. [I143II and H 14411 that the "cut-off state" of the corresponding coherent 
state, 

«Aac,/A) = .a;c,/A)/||<a;c,/A)||, (97) 
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with 



s,y(s)=y 

has desired semiclassical properties in testing the kinematical operators (e.g. holonomy and 
flux). But unfortunately, these cut-off coherent states cannot be directly used to test the semi- 
classical limit of the Hamiltonian constraint operator S{N), since S{N) is graph-changing so that 
its expectation values with respect to these cut-off states are always zero! So further work in this 
approach is expected in order to overcome the difficulty. Anyway, the complexifier approach 
provides a clean construction mechanism and manageable calculation method for semiclassical 
analysis in loop quantum gravity. 

We now turn to the second approach. As we have seen, loop quantum gravity is based on 
quantum geometry, where the fundamental excitations are one-dimensional polymer-like. On 
the other hand, low energy physics is based on quantum field theories which are constructed 
in a flat spacetime continuum. The fundamental excitations of these fields are 3-dimensional, 
typically representing wavy undulations on the background Minkowskian geometry. The core 
strategy in this approach is then to relate the polymer excitations of quantum geometry to Fock 
states used in low energy physics and to locate Minkowski Fock states in the background inde- 
pendent framework. Since the quantum Maxwell field can be constructed in both Fock repre- 
sentation and polymer-like representation, one first gains insights from the comparison between 
the two representations, then generalizes the method to quantum geometry. A "Laplacian oper- 
ator" can be defined on 'Hkin [.26.1 1,191. from which one may define a candidate coherent state Oq, 
also in Cyl* , corresponding to the Minkowski spacetime. To calculate the expectation values of 
kinematical operators, one considers the so-called "shadow state" of Oq, which is the restriction 
of Oo to a given finite graph. However, the construction of shadow states is subtly different from 
that of cut-off states. 

We will only describe the simple case of the Maxwell field to illustrate the ideas of the 
construction [|152||[|153ll l|20l. Following the quantum geometry strategy discussed in Sec. 4, the 
quantum configuration space A for the polymer representation of the U{\) gauge theory can be 
similarly constructed. A generalized connection A 6 A assigns each oriented analytic edge in 
2 an element of U{\). The space A carries a diffeomorphism and gauge invariant measure //q 
induced by the fiaar measure on U{1), which gives rise to the fiilbert space, 'Ho := L^(A, diJ.o), 
of polymer states. The basic operators are holonomy operators A(e) labeled by one-dimensional 
edges e, which act on cylindrical functions by multiplication, and smeared electric field opera- 
tors E{g) for suitable test one-forms g on S, which are self-adjoint. Note that, since the gauge 
group U{\) is Abelian, it is more convenient to smear the electric fields in 3 dimensions [|20ll . 
The eigenstates of E{g), so-called flux network states Na,n, provide an orthonormal basis in 'Hq, 
which are defined for any finite graph a with edges as: 

Na,n{^) := [A(ei)]"' [A(e2)]"^ • • • [A(e;v)]"^ (99) 

where ft = (ni, ■ ■ s^a^) assigns an integer ni to each edge e/. The action of E{g) on the flux 
network states reads 

E{g)N,^n = -KYj^i r g)Nafl. (100) 
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In this polymer-like representation, cylindrical functions are the finite linear combinations of 
flux network states and span a dense subspace of "Hq. Denote by Cyl the set of cylindrical 
functions and by Cyl* its algebraic dual. One then has a triplet Cyl c "Ho c Cyl* in analogy 
with the case of loop quantum gravity. 

The Schrodinger or Fock representation of the Maxwell field, on the other hand, depends on 
the Minkowski background metric. Here the Hilbert space is given by "Hp = L^{S' , dup), where 
S' is the appropriate space of tempered distributions on E and is the Gaussian measure. The 
basic operators are connections A(/) smeared in 3 dimensions with suitable vector densities 
/ and smeared electric fields E{g). But A(e) fail to be well defined. To resolve this tension 
between the two representations, one proceeds as follows. Let x be the Cartesian coordinates of 
a point in 2 = R^. Introduce a test function by using the Euclidean background metric on R^, 

/^(^ = n 111 3 exp(-|:^V2r2), (101) 

which approximates the Dirac delta function for small r. The Gaussian smeared form factor for 
an edge e is defined as 

X° ,)(f) := r dsms) - (102) 
Then one can define a smeared holonomy for e by 

A(,)(e) := exp[-/ r X^,,)(f)A„(f)], (103) 
Jr3 

where Aa{x) is the ?7(1) connection one-form of the Maxwell field on 2. Similarly one can 
define Gaussian smeared electric fields by 

£wte):= [ gaix) r Sriy - ^E'^iS). (104) 
Jr3 Jr3 

In this way one obtains two Poission bracket algebras. One is formed by smeared holonomies 
and electric fields with 

{A(,)(e),A(,)(e')} = = {E{,g\E{g')\ (105) 

{A(,)(e),£(g)} = -/( r Xl^^ga)Kr)ie). 
Jr3 

The other is formed by unsmeared holonomies and Gaussian smeared electric fields with 

{A(e), A(e')} = = {E(,)(g), E(,)(g')} (106) 
{A(e),£(,)(g)} = -/( r X" ,)g,)A(e). 

JR3 

Obviously, there is an isomorphism between them, 

/, : (A(,)(e), E{g)) ^ (A(e), E^^yig)). (107) 
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Using the isomorphism /,, one can pass back and forth between the polymer and the Fock 
representations. Specifically, the image of the Fock vacuum can be shown to be the following 
element of Cyl* [fT52ll[fT53ll . 



iV\ = J] exp[-- ^ Gjjmnj] (NU (108) 

a,ft IJ 

where {Na,n\ £ Cyl* maps the flux network function \Na,it) to one and every other flux net- 
work functions to zero. While the states {Na,n\ do not have any knowledge of the underlying 
Minkowskian geometry, this information is coded in the matrix Gu associated with the edges 
of the graph a, given by GOll 

G,j = f dte%t) r dt'ej\t') r d'x6ab{^Vfri.x-Ut)mVh{xJj{t'))l (109) 

Jei Jej J 

where 6ab is the flat Euclidean metric and A its Laplacian. Therefore, one can single out the Fock 
vacuum state directly in the polymer representation by invoking Poincare invariance without 
any reference to the Fock space. Similarly, one can directly locate in Cyl* all coherent states 
as the eigenstates of the exponentiated annihilation operators. Since Cyl* does not have an 
inner product, one uses the notion of shadow states to do semiclassical analysis in the polymer 
representation. From Eg. (11081) . the action of the Fock vacuum {V\ on Na,n reads 

(V\Na,H) = r dixlV,Na,n, (110) 
Ja„ 

where the state Va is in the Hilbert space "Kq, for the graph a and given by 

y,(A) = e^P[-? Yj Na,n{P^). (HI) 

Thus for any cylindrical functions ipa associated with a, 

{V\ip,) ={VM, (112) 

where the inner product in the right hand side is taken in "K^. Hence Va{^) are referred to as 
"shadows" of {V\ on the graphs a. The set of all shadows captures the full information in {V\. 
By analyzing shadows on sufficiently refined graphs, one can introduce criteria to test if a given 
element of Cyl* represents a semi-classical state [|20ll . It turns out that the state {V\ does satisfy 
this criterion and hence can be regarded as semi-classical in the polymer representation. 

The mathematical and conceptual tools gained from simple models like the Maxwell fields 
are currently being used to construct semiclassical states of quantum geometry. A candidate 
kinematical coherent state corresponding to the Minkowski spacetime has been proposed by 
Ashtekar and Lewandowki in the light of a "Laplacian operator" [|T9]| [1201 . However, the detailed 
structure of this coherent state is yet to be analyzed and there is no a priori guarantee that it is 
indeed a semiclassical state. 

One may find comparisons of the two approaches from both sides H 14711 ll20l . It turns out that 
Varadarajan's Laplacian coherent state for the polymer Maxwell field can also be derived from 
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Thiemann's complexifier method. However, one cannot find a complexifier to get the coherent 
state proposed by Ashtekar et al. for loop quantum gravity. Both approaches have their own 
virtues and need further developments. The complexifier approach provides a clear construc- 
tion mechanism and manageable calculation method, while the Laplacian operator approach is 
related closely with the well-known Fock vacuum state. One may also expect that a judicious 
combination of the two approaches may lead to significant progress in the semiclassical analysis 
of loop quantum gravity. 

6.2 Algebraic Quantum Gravity Approach 

As we have shown in the last subsection, although Thiemann's complexifier coherent state has 
a clear calculable mechanism and correct semi-classical properties in testing kinematical opera- 
tors, it fails to be a qualified semi-classical state for the quantum dynamics since the semiclassi- 
cal limit of the Hamiltonian constraint operator S{N) or master constraint operator M is clearly 
not correct, both S(N) and M are graph-changing so that their expectation values with respect 
to these cut-off coherent states are always zero. So a possible way to avoid such a problem is to 
define a non-graph-changing version of Hamiltonian constraint operator or similarly, a master 
constraint operator. However, such a modification is hard to make in the framework of loop 
quantum gravity since the action of Hamiltonian constraint operator always adds several arcs 
on certain graphs. But if the framework of loop quantum gravity is suitably modified then it 
turns out that a version of non-graph-changing Hamiltonian constraint operator can be proposed 
and the semi-classical analysis for the quantum dynamics can be carried out with the complex- 
ifier coherent states defined previously. Such a modification is recently made by Thiemann in 
ll79ll[|80l[|8T| and is called algebraic quantum gravity (AQG) approach. We describe it briefly in 
what follows. 

Algebraic quantum gravity is a new approach to canonical quantum gravity suggested by 
loop quantum gravity. But in contrast to loop quantum gravity, the quantum kinematics of al- 
gebraic quantum gravity is determined by an abstract *-algebra generated by a countable set 
of elementary operators labeled by a single algebraic graph with countably infinite number of 
edges, while in loop quantum gravity the elementary operators are labeled by a collection of 
embedded graphs with finite number of edges. Thus one can expect that in algebraic quantum 
gravity, we lose the information of the topological and differential structure of the manifold in 
all the quantization procedure before we do semi-classical analysis. Hence the quantum the- 
ory will be of course independent of the topology and differential structure of the manifold but 
based only on an algebraic graph, which only contains the information of the number of vertices 
and their oriented valence. 

Definition 6.2.1: An oriented algebraic graph is an abstract graph specified by its adjacency 
matrix a, which is an N x N matrix. One of its entries ajj stand for the number of edges that 
start at vertex I and end at vertex J. The valence of the vertex I is given by V/ = Yjji^^u + ^ 
We also use V{a) and E{a) to denote the sets of vertices and edges respectively. 

In our quantization procedure, we fix a specific cubic algebraic graph with a countably infi- 
nite number of edges N = ^ and the valence of each vertex v/ = 2 x dim(E). Such a specific 
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choice, although it detracts from the generality of the theory, is practically sufficient for our use 
in the semiclassical analysis. 

Given the algebraic graph a, we define a quantum *-algebra by associating with each edge 
e an element A{e) of a compact, connected, semisimple Lie group G and an element Ej{e) take 
value in its Lie algebra Q. These elements are subject to the commutation relations 

[A(e),Aie')] = 0, 

[Ej{e),A(e')] = inQ\yTj/2A{e), 
[Ej(e),Aie')] = -ifiQ^6e,e'fjkiEi{e'\ 

and * -relations 

A{er = {A{er'Y, Ejie)* = Ej(e), 

where Q stands for the coupling constant, tj is the generators in the Lie algebra g and fji^i is the 
structure constant of g. We denote the abstract quantum *-algebra generated by above elements 
and relations by %. 

A natural representation of % is the infinite tensor product Hilbert space "K® = ®e'7Y^ where 
He = L}(G, dfifj) H45h whose element is denoted by ®f = ®efe- Two elements ®f and ®f> in 
"H® are said to be strongly equivalent if 2e I < fe, f'e >'H, -1| converges. We denote by [/] the 
strongly equivalence class containing (8)/. It turns out that two elements in "K^ are orthogonal 
if they lie in different strongly equivalence classes. Hence the infinite tensor Hilbert space "K^ 
can be decomposed as a direct sum of the Hilbert subspaces (sectors) "K^^j which are the closure 
of strongly equivalence classes [/]. Furthermore, although each sector "K^^j is separable and 
has a natural Fock space structure, the whole Hilbert space 'H^ is non-separable since there are 
uncountably infinite number of strongly equivalence classes in it. Our basic elements in the 
quantum algebra are represented on 'H^ in an obvious way 

A{e)®f := [A(e)fA®[®e'tefe'], 
Ej(e)®f := [inQ^X'jfe]®[®,>^Je'l 

As one might have expected, all these operators are densely defined and Ej{e) is essentially 
self-adjoint. Given a vertex v 6 V{a), the volume operator can be constructed by using the 
operators we just defined 

K := £l Jl^ Yj ^"(^1,^2, e3)e'J'^Ei{ei)Ej(e2)Me3)\, 

where the values of €^.{61,62,63) should be assigned once for all for each vertex. When we 
embed the algebraic graph into some manifold, the embedding should be consistent with the 
assigned values of 6^(61,62, e^). 

Then we discuss the quantum dynamics. By the regularization methods frequently used in 
the last two sections, the half densitized constraints can be quantized to be composite operators 
as we list below. 
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• Gauss constraint 

Gj(v) := <2t'/'> J] Ejie); 

e at V 

• Spatial difFeomorphism constraint 

^ E(v) ^ \L(v,eue2)\ 



Euclidean Hamiltonian constraint (up to an overall factor) 



6,(g|,^^2,g3) 

|L(v,ei,e2)l 



X 2 Tr([AOe) - A{fir']A{e^){A{e,r\ V^'^])- 

jS€L(v,ei,e2) 

• Lorentzian Hamiltonian constraint (up to an overall factor) 

•= TTT y! 6v(ei,e2,e3) 
£(v) 

X Tr(A(ei)[A(ei)-\ [^^^\ y]]A(e2)[A(e2)-\ {A{e^)-\ {Uf, V]] 

X A(e3)[A(e3)-\ V^]), 
H{v) = H^^l'-\v) + f{v)- (113) 

where V := K, := Zv^^^W and 

•= X 6v(ei,e2,e3) 

£^(v) 

X Tr(A(ei)[A(ei)-i, yW]A(e2)[A(e2)"\ t/^'■^]A(^^3)[A(^^3)"^ V^'')]). 

L(v, ei , ^2) denotes the set of minimal loops starting at v along ci and ending at v along e~^. And 
a loop yS e L(v, ^1,62) is said to be minimal provided that there is no other loop within a sat- 
isfying the same restrictions with fewer edges traversed. Note that since we only have a single 
cubic algebraic graph, the difFeomorphism constraint can only be implemented by defining the 
operators corresponding to diff'eomorphism generators because a finite difFeomorphism trans- 
formation is not meaningful in our algebraic treatment unless the algebraic graph is embedded 
in a manifold. As a result, the (extended) master constraint can be expressed as a quadratic 
combination: 

v€V(a) 

76 



It is trivial to see that all the above operators are non-graph-changing and embedding indepen- 
dent because we have only worked on a single algebraic graph so far. However, when we test 
the semiclassical limit of these operators, especially the master constraint operators, we should 
specify an embedding map X which map a algebraic graph to be an embedded one. With this 
specific embedding, we can see the correspondence between the classical algebra of elementary 
observables and the quantum *-algebra. We define the holonomy and suitably modified flux by 

A(e) := AiXie)) :=!Pexp( f A), 

JX{e) 

Ej{e) := -2Tr[Ty f eatcdx" A dx'A{p,ix))E'ix)A{p,(x))-'l 
Js, 

where 5 ^ is a face which intersects the edge X(e) only at an interior point pe of both S e and X(e). 
We choose a system of paths {peix)}^ for all x e S e, such that Peix) starts at s(X(e)) along X(e) 
until Pe and then runs within 5 ^ until x. As one might expect, the quantum *-algebra we defined 
previously is just consistent with the classical Poisson algebra generated by these holonomis 
and fluxes: 

{A(e),A(e')} = 0, 
{Ej{e),A{e')} = Q^SeyTj/2A(e), 
{Ej{e\A{e')} = -Q'5e,e'fjkiEi{e'). 

Then we consider the coherent states. By employing the Laplacian complexifier on each 
edge 

1 

the coherent state is obtained as it was in the last section: 



n 

where A„ denotes the eigenvalue of the Laplacian on G and t = ip/aj represents the classical- 
ization parameter. The coherent state peaks at the complexified classical phase space point 

g(A,E) := y ^-J-{Ce,A(e)], = exp(/£(e)/a^)A(e), 

note that the parameter Ue is specified such that E{e)/al is dimensionless. Hence the coherent 
state on the whole graph is represented by an infinite tensor product state: 



eeE{a) ' ' e;(A,E) 



The peakness, fluctuation and other semiclassical properties of these states have been checked 
in [I143II HI 441 in which the most important part is that 

< ^A,E\A(eWA,E >= A(e) < ¥A,£|£(e)|YA,£ >= E{e) 
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up to terms which vanish faster than any power of te as te 0. And the fluctuations are small. 

With the semiclassical state we just constructed, the expectation value of the above (ex- 
tended) master constraint operator can be calculated and its semiclassical limit can be tested. In 
the following, we summarize the result of the calculation. In [801, a semiclassical calculation 
for the master constraint operator is carried out based on a cubic algebraic graph. The calcula- 
tion makes use of a simplifying assumption: we substitute the gauge group for gravity S U{2) 
by U(l)^. And the result of the calculation shows that in U{1)^ case the (extended) master 
constraint operator has correct semiclassical limit 

lim < XM\X« >= M^^^^Lm] ^ M[m] (e 0) 

where m represents a phase space point and e is the lattice parameter such that the lattice be- 
come continuum as e ^ 0. In addition, it is shown that the next-to-leading order terms which 
contribute to the fluctuation of M are finite. 

Moreover, the calculation in [81 J shows that the result of the exact non-Abelian calculation 
matches precisely the results of the Abelian approximation, provided that we replace the clas- 
sical f/(l)^ terms {H, p'j}j=ix3 by {Tr(TjA(e)),Tr(TjE(e))}j=ix3, which means that the theory 
of algebraic quantum gravity admits a semiclassical limit whose infinitesimal gauge symmetry 
agrees with that of general relativity. 
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7 Conclusion and Discussion 



As it was shown in the previous sections, loop quantum gravity offers a conceptually clear 
and mathematically rigorous approach to quantize general relativity. In this approach, we are 
seeking new physics deeply below the Planck scale. In the kinematical framework, a quantum 
Riemannian geometry is established and some geometrical operators, e.g. area, volume, are 
well-defined, and their spectrum are shown to be discrete, which means that the structure of the 
space may be discrete below the Planck scale. Such a new phenomena sheds light on quantum 
field theory, lattice gauge theory and their renormalization. Moreover, the program in the quan- 
tum dynamics of loop quantum gravity represents significant progress in the research area of 
quantum gravity. Before loop quantum gravity, the quantum Wheeler-DeWitt equation was only 
a formal equation and from concrete calculations. However, in the framework of loop quantum 
gravity, we already have a well-defined quantum Hamiltonian constraint operator which has 
an explicit action on kinematical states, so that the quantum Wheeler-DeWitt equation is well- 
defined in loop quantum gravity. On the other hand, the matter field can also be quantized in 
this framework and we show that the matter Hamiltonian is free of UV-divergence and don't 
need a renormalization process. Furthermore, with the coupled matter field, a matter coupled 
Hamiltonian constraint operator is obtained so that the problem of time may be solved and, such 
an idea is being translated into a new understanding of the early universe in the context of loop 
quantum cosmology. 

Although great progress has been made, as an unfinished framework, loop quantum gravity 
still has many issues to be solved in the future research. To conclude this thesis, we list some of 
those in the following: 

• First of all, we don't have the complete solutions for either Hamiltonian constraint equa- 
tion or master constraint equation. Thus one cannot explicitly construct the physical 
Hilbert space for loop quantum gravity. So the quantum dynamics of gravity is essen- 
tially unknown so far. 

• To make contact with experimental results, one should know the observables in the quan- 
tum theory which have to be invariant under gauge transformation. However, some of the 
Dirac observables that have been constructed involve an infinite number of derivatives 
and extremely hard to manage [il50ll[[6T1|[|62|| . 

• The semiclassical limit of loop quantum gravity is unknown so far, although a great deal 
of progress has been made in the context of algebraic quantum gravity. And in alge- 
braic quantum gravity, further research work is needed to show the fluctuation of master 
constraint operator should be small. 

• As it was shown at the end of section 4.3, the regularization process for the Hamiltonian 
constraint operator is ambiguous and there is a list of free parameters. Thus it is also a 
research project to remove as many ambiguities as possible. And some work has been 
recently done in this direction UllOL 

• The Immirzi parameter is another free parameter in the framework of loop quantum grav- 
ity, which comes in with the classical formulation. In the classical theory, diff'erent values 
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of the Immirzi parameter label equivalent classical theories since they are connected by 
canonical transformations. However, in quantum theory, it is problematic because the 
representation with different Immirzi parameter are not unitarily equivalent. 

The construction of loop quantum gravity crucially depends on the compactness of its 
gauge group SU(2), which comes from an internal partial gauge fixing. And it is ar- 
gued that the internal Lorentz symmetry is broken in a non-natural way [I131II . So it 
seems to be better to switch back to the complex Ashtekar variables which are free 
of the internal gauge fixing and preserve the internal Lorentz symmetry and, will also 
greatly simplify the Hamiltonian constraint. However, the price is that we should work 
on non-compact gauge group 5L(2, C), and there is no satisfactory quantization pro- 
gramme for the 5L(2, C)-gravity so far, although some work has been done in this di- 
rection [l7T]|[fT05ll[fro6ll . 

As it was shown in section 6.2, algebraic quantum gravity provides a clear way to make 
many calculations accessible and the semiclassical analysis can be carried out in this 
framework. And it seems that the construction of algebraic quantum gravity admits the 
non-compact internal gauge group (or non-compact reduced configuration space) since 
there is only one graph in quantization process. However, in this framework, the restric- 
tion of the possible representation is so loose that there even exists a possible representa- 
tion in which the spectrum of geometrical operators are continuous. 

The transition amplitude calculation for loop quantum gravity is accessible in the so- 
called spin foam model and depends on a unclear conjecture of GFT/spinfoam duality. 
On the other hand, it is still not clear how to build a path-integral formulation and connect 
with spin foam models from the canonical approach, and such a connection may give the 
needed support for GFT/spinfoam duality. 

We have constructed the dynamics of matter quantum field theory on a quantum back- 
ground in section 5. However, it is not clear how we can make the connection with the 
ordinary quantum field theory on curved spacetime. Moreover, it is still an problem how 
to find the non-perturbative correspondence of Hadamard states in perturbative quantum 
field theory in curved spacetime, although some hints of a connection may come out in 
spin foam calculations [|127ll[|42ll . 
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